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FOREWORD 


The research described in this report was performed for the 
Laser Systems Section at the NASA/Lewis Research Center 
under Contract NAS3-20381. The NASA Project Manager was Mr. 
Stephen M. Cohen. 

The beam shape used in Subsection 4. 2 was suggested by Dr. 
Peter Nebolsine. Efficient computational assistance was provided by 
Edward Tomaszewski. 


SUMMARY 
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Reference 1 describes a study of a CW laser-heated hydrogen rocket 
using 10. 6 ^rn radiation. It was found that the properties of Laser Supported 
Combustion (LSC) waves in hydrogen were an important ingredient in model- 
ing this rocket. They determine the temperature level reached by the gas, 
and the mass flux through the rocket. No theoretical or experimental studies 
of LSC waves in hydrogen were available. Therefore, a theoretical study of 
such waves has been made. A one -dimensional energy equation, with con- 
stant pressure and area, was used to model the LSC wave. This equation 
balances convection, conduction, laser energy absorption, radiation energy 
loss and radiation energy transport. The latter is shown to be approximated 
well by a radiation conduction model. Solutions of this energy equation were 
obtained to give profiles of temperature and other properties, as well as 
the relation between laser intensity and mass flux through the wave. The 

calculations cover the pressures of 1, 3, 10 and 30 atm, 10. 6 |jm 

4 2 6 2 

laser intensities from 10 W/cm to 10 W/cm and power levels of 10 kW 
and 5 MW. The physics of these waves leads to high peak temperatures (on the 
order of 20, 000K), because the absorption mechanism is inverse bremsstrahlung, 
which requires a significant degree of ionization, and this occurs in hydrogen 
only above about 10, 000K. The high temperatures also lead to considerable 
radiation losses. 

The flow through the LSC wave was then conducted through a variable- 
pressure, variable-area streamtube to accelerate it to high speed, with the pro- 
pulsion application in mind, A numerical method for coupling the LSC wave 
model to the streamtube flow was developed, and a sample calculation was 
performed. The result shows that 42% of the laser power has been radiated away by 
the time the gas reaches the throat. This is in contrast with results for a 
similar case calculated in Ref. 1, where the LSC wave properties were only 
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estimated, not calculated. There only 5% of the power was found to be lost. 
The present more realistic calculations show the large losses incurred by 
the necessity of operating at high temperatures in hydrogen. 

Several two-dimensional effects were estimated. Radial losses due 
to heat conduction and black radiation from the edge of the hot gas were found 
important at the 10 kW power level. They cause the intensity threshold for 
existence of the LSC wave to rise to values considerably higher than that 
which prevails if these losses are ignored. Beam convergence effects were 
also estimated by including beam area change in the laser absorption term. A 
calculation showed the mass flux required in a converging beam to be larger 
than in a parallel beam of the same initial intensity. Converging beams pro- 
vide static stability for the wave position. Transverse velocity was also con- 
sidered. It was concluded that in the radially confined flows of interest for 
propulsion applications, transverse velocities would be less important than 
in the unconfined flows where air experiments have been conducted. 

It would be advantageous to produce LSC waves at lower temperatures 
to reduce the large radiation losses. This can be accomplished by introduc- 
ing easily ionized seed into the hydrogen to allow laser absorption at lower 
temperatures. The amount of seed must be small enough so that the specific 
impulse is not adversely affected by its weight. It is recommended that 
studies of LSC waves in seeded hydrogen be conducted. 
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1. INTRODUCTION 


In a recent study of CW laser propulsion. Ref. 1, an analysis 
was made of the specific impulse which could be obtained by heating 
hydrogen with a 10. 6 ^m laser beam. With this combination of gas 
and laser frequency, the only mechanism by which the laser energy 
can be absorbed ; n the gas is inverse bremsstrahlung. This mechanism 
requires the presence of electrons and therefore will only become 
operable when the hydrogen ionizes at temperatures of 10, 000K or 
more. Other phenomena must operate to heat the gas from ambient tern- 
perature to these high temperatures. 

It has been known for some time that a mechanism to achieve this 
heating exists for air, known as a laser- supported combustion (LSC) wave. 
This wave is analogous to a deflagration wave in ordinary chemical com- 
bustion, but the energy source is a laser beam, rather than chemical 
energy. The LSC wave consists of an approximately constant pressure 
region in which the gas temperature rises sharply to a high value, and 
then decays slowly. The high temperature region is maintained by the 
absorption of the laser energy by the gas. The rapid temperature rise 
is produced by the transport of energy from the hot region into the cold 
gas by both thermal conduction and radiation. The decay of the temp- 
erature at the back of the wave is caused by radiation losses dominating 
the absorption of laser energy, because the energy in the laser beam 
has already been almost completely absorbed by the gas in the hot 
region. 

The energy transport by conduction and radiation in an LSC wave 
therefore provides a mechanism to heat cold gas up to a temperature where 
laser energy can be absorbed by a high temperature mechanism like inverse 
bremsstrahlung. Such waves have been observed in air and the rare 
gases, though there are no reports of any attempt to generate them in 
hydrogen. 
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The laser-heated hydrogen rocket is thus envisioned to have an 
LSC wave standing at the front of the heating region. Cold hydrogen 
flows in, is heated and absorbs laser energy at high temperatures. It 
is then expanded through a nozzle to a high velocity to impart thrust to 
the rocket. 

Raizer (Ref. 2) was the first to recognize the possibility of 
such a wave, though he considered primarily thermal conduction, not 
radiation, as the forward energy transport mechanism. 

Boni and Su (Ref. 3) and Jackson and Nielsen (Ref. 4) added the 
effect of radiation, which they found to be important in many cases. 
Additional modeling for air has been done Boni and Su (Refs. 5 and 6) 
and by others. Recently Pirri et al (Ref. 7) have improved the radiation 
model in an air LSC wave, and developed a simple numerical calculation 
scheme for computing the temperature profile. 

Experimental observations of LSC waves in air have been made 
by Klosterman and Byron (Ref. 8) among others, and in rare gases by 
Soviet investigators, but not in hydrogen. 

An important feature of these waves is their speed. The models 
indicate that the speed of an LSC wave in a given gas is a unique function 
of the laser intensity. Since the wave in the laser-heated rocket is 
to be stationary, it is important to know the relationship between speed 
and intensity, so the proper inflow velocity of the hydrogen to hold 
the wave fixed may be chosen for the intensity used. 

This inflow determines the mass flow per unit area in the 
rocket, and so affects the thrust level which can be achieved for a given 
cross- sectional area. Therefore, thrust and power scale linearly with area. 

In the study of Ref. 1 it was noted that in the absence of either 
models or experiments for hydrogen, an estimate would have to be 
made of wave speed versus intensity. Such an estimate was made, 
using the method of Ref. 4 and a crude approximation to the radiation 


- 2 - 


properties of hydrogen. However, the uncertainty in the resulting 
speed was felt to be considerable. Since this parameter is so important 
to the design of a laser-heated hydrogen rocket, and since experiments 
on LSC waves and rockets are underway at NASA Lewis Research 
Center, it seems desirable to develop a model of LSC waves in hydrogen. 
The work reported here describes such a model, based on the com- 
paratively simple numerical methods of Ref. 7, and detailed considera- 
tion of the thermodynamic and radiation properties of hydrogen. 

In Section 2, a one-dimensional model of the hydrogen LSC 
wave is developed by solving the energy conservation equation, including 
convection, laser absorption, radiation losses to the environment, 
radiation transport of energy from one part of the wave to another, and 
thermal conduction. Calculations are made for laser Intensities from 
10^ to 10^ W/cm^, at pressures of 1, 3, 10 and 30 atm, for laser powers 
of 10 kW and 5 MW. The resulting wave speed, mass flow and temp- 
erature parameters are tabulated for 45 cases. 

In Section 3, the acceleration of the gas coming out of the LSC 
wave by a variable area, variable pressure channel is considered. 

This represents the nozzle necessary to accelerate the gas to produce 
thrust. The purpose of this section is to develop a computational method 
for coupling the LSC wave model to the nozzle flow. One case is com- 
puted to demonstrate the use of the proposed method of coupling for a 
specific set of parameters. 

Section 4 is concerned with two-dimensional effects which may 
be of importance to the wave speed. We consider radial heat losses 
due to conduction and black radiation emission from the edges of the hot 
gas cylinder; convergence of the laser beam; and transverse flow 
velocities. 

Section 5 sets forth the conclusions which can be drawn, and 
how they affect the results of Ref. 1. 
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There are also four Appendices, which contain details of 
various aspects of the model. Appendix A gives the model for the 
thermodynamic properties of equilibrium hydrogen used. This repre- 
sents an improvement over the similar model described in Ref. 1, 
because it includes the lowering of the ionization potential, and a 
variable cut-off for the electronic partition function of the hydrogen 
atom. 

Appendix B defines the absorption coefficient for 10.6 lim 
radiation in ionized hydrogen by inverse bremsstrahlung. 

Appendix C gives the values for the thermal conductivity of 
equilibrium hydrogen used. 

Appendix D describes briefly th radiation properties of hydrogen, 
and gives the expressions used for axial radiation transport of energy 
and for loso of energy by radiation. 
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2. ONE DIMENSIONAL LSC WAVE MODEL 


2. 1 Basic Equations 

The LSC wave is described by an energy balance for the flowing gas 
and the incoming laser energy. The elements that make up this balance 
include absorption of laser energy, thermal conduction of energy, transport 
of energy by radiation, convection of energy, and loss of energy by radia- 
tion. The equation which expresses ’■bis balance is 


_d 

dx 


j» d S. 

. TIV . dl A 

(puH) t— t — 


+ V£ 



( 2 . 1 ) 


Here p is gas density, u is gas velocity, H is gas total enthalpy (h + u /2) 

where h is the thermochemical enthalpy), I is the laser intensity, 

S the radiation energy flux, P the radiation loss term, and X the equili< 
AT ^ 

brium thermal conductivity of the gas. The terms in this energy equation 
are, respectively, convection, laser energy absorption, radiation energy 
transport, radiation energy loss and thermal conduction. 

Conservation of mass is satisfied by the statement that the mass 
flow per unit area is a constant, 


pu ’ a a 


= constant 


( 2 . 2 ) 


Conservation of momentum is replaced by a constant pressure 
condition. 


p = Po 


(2.3) 
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because the velocity changes are so small that pressure changes generated 

by momentum changes are negligible. 

The laser absorption term is obtained from the laser energy trans- 
port term 


dl/dx = - k I 

w h e re k is the absorption coefficient for the laser energy. The solution 

l_j 

of this equation is 


= I o ;T ' T= / k L 


dx 


(2.4) 


where 


T is the optical depth and I q is the incoming laser flux. Thus we find 


dl/dx = - k L I o e 


( 2 . 5 ) 


The convection term, using (2. 2), becomes 


d . /dh, du\ 

^(puH) = m A ^ + u-) 


. ( dT + ^d £ " 1 \ 

" m A yp dx p dx / 


( 2 . 6 ) 


The equation of state, (A-8), provides p in terms of T for the constant 


pressure p 
r o 


-1 


P = R M ZT/P o 


(2.7) 
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where R M is the molecular gas constant, and Z is the compressibility factor. 
Differentiation gives 


5aa 








•?/,* a'* • ,’^r* * 


-(PUK, = & 


Z 2 T L t 41 

pl V dfnT / dx 


( 2 . 8 ) 


As discussed in Appendix D, we will use the conduction approximation 
for the radiation transport so that 


A = _ _d 
dx dx 


(vS) 


(2.9) 


where X is radiation conductivity. 
K 


Use of (2.5 ), (2. 8) and (2. 9) enables us to express the energy equation 




A -SS S 


1 r-T.-r, d m , , . dT 

“ k L I o® P T dx (X C + X R ) dx 


( 2 . 10 ) 


The quantities Z, dinZ/d^nT, c , k , X , X n and P are given in the 

p J-j K 1 


Appendices as functions of T for a given p = p . They appear in Eqs. (A- 19) 

o 


(A-20), (B-4), (C-i) and Table Cl, (D-9) and (D-l). Then (2. 10) is a 
second order differential equation for T(x), the temperature distribution in 
the LSC wave. When suitable boundary conditions are imposed, it can be 


solved. 
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One set of boundary conditions can be obtained from the incoming 
flow at x = - ®. However, one cannot solve the equations numerically start- 
ing at the incoming conditions, because these conditions are a solution to the 
equation and the numerical procedure will not move off that solution. Instead, 
one can integrate (2. 1) from x = - ® where the initial conditions prevail, to 
x = 0, which is an arbitrary location, except that our definition of t, (2.4), 
provides that no absorption of laser energy has yet occurred there. Denoting 
conditions at x = 0 by subscript i, we have 


m [h. - h + (u 2 - u 2 )/2] 
A 1 i o i o J 


0 

+ / P T d *= ft C + V 
— 00 


/ dT \ 

i Vdx/. 
' ' l 


The term on the right side at x = - • has been set to zero since the tempera- 
ture gradient is zero there. 

We now assume that the temperature T. is low enough at x = 0 so 
that no radiation losses have occurred yet, thus ignoring the integral over 
P p. Replacing u by the mass conservation and state equations yields 


. T . 7* a r m\ 

m. h. - h +1 I 

A L i o \ p / 


(Zf T - T ) 
11 o 


] ■ ftc + Vi (fa) 
J ' '1 


( 2 . 11 ) 


The enthalpy in (2. 11) is given as a function of T in (A- 27). For a chosen 
initial temperature T., (2. 11) specifies (dT/dx)., and these give the two 
boundary conditions necessary to integrate the second order equation (2. 10). 
The result should be independent of the choice of T., as long as it is low 
enough so that no significant laser energy absorption and no significant 
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radiation energy loss have yet occurred* This can be checked by performing 
the integrations for several values of T. and verifying that the results do 
not change. 


2. 2 Scheme for Numberical Solution 

The integration requires the specification of the mass flow constant 

m for a given laser flux I . As previous studies of LSC waves have shown, 
A ° 

the integration of the energy equation is an eigenvalue problem. There is 
only one value of m A for a given I which allows a physically realistic 
solution. This is related to the downstream conditions at the back of the 
LSC wave, as shown in Fig. 2. 1. If m^ is too high, say m^^, the tempera- 
ture profile will either rise indefinitely or peak, fall, and then rise again. 
This indicates that there is so much mass flow that the laser energy cannot 


heat it all, and an energy source is present downstream. On the other hand, 
if the mass flow is too low, say ™ AL , the temperature profile will rise to 
a peak, and then fall rapidly to zero, indicating the presence of an energy 
sink downstream. Only one value of m^ will permit the temperature to 
reach a finite value far downstream, and the integration must be performed 
so as to find this value* 

We use an iterative scheme for this purpose* First, integrations are 
performed for several values of m f until a pair of bounds are found, of 
which the upper bound causes the temperature to rise again after peaking and 
decreasing, and the lower bound causes the temperature to fall continuously 
after the peak. These two are then used in a bounding scheme. Their aver- 
age is used for an integration. If it behaves as m^y, it is taken as the new 
upper bound, and if it behaves as m AL it becomes the new lower bound. 

We continue to perform integrations, always using the average of the closest 
upper and lower bounds as the next value of Eventually, an( * 

agree to as many significant figures as we like, and the corresponding 
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temperature profiles agree very closely to some value of x, after which 

they diverage as shown in Fig. 2. 1. The common value of m * m 

. AU AL 

and the common temperature profile is the solution sought. 

Notice that we never are able to attain a profile which approaches 
a finite asymptotic temperature as x -* ». This would require infinite pre- 
cision in m A> However, in our model there is no real significance to 
this asymptote. If we could find it, its value would depend on the low- 
temperature details of the radiation energy transport and loss, which are 
not accounted for correctly in our model anyway. Our procedure provides 
m A' tlie tem P erature profile in the region of high temperature, the absorp- 
tion of laser energy, and the radiation losses, which are the important 
features of the LSC wave. The low temperature tail has no real significance. 

The integration is performed using a standard 4th order Runge- 
Kutta integration scheme for a system of first order equations. The present 
system of equations is written in that way by defining a new variable W by 

dT _ W 

d x- V X R (2.12a) 


Then (2. 10) is 


dW _ m A W 
dx X , 


Sv Ft?) «■■(■•») 


- k I e T + P 
L o T 


(2. 12b) 
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where, from (2.4) 


dT 

dx 



(2. 12c) 


This is a third order system for T, W, T, whose boundary conditions are 

x = Os T = T\, T ~ 0 (2. 13a) 


and from (2. 14) 



2 2 
T 

l i 



(2. 13b) 


The parameters necessary to integrate this system are I q » m^, I\, P q , T q , 

and the radius R . of the plasma tube. This latter parameter enters because 

the radiation loss P T and the radiation conduction coefficient X R depend on 

R . as described in Appendix D. R. is found by choosing the laser power P 
1 * 
and relating it to the intensity and area of the plasma tube as 


P = I tt R 2 , R. = (P/I IT) 1/2 
o i i ° 


(2. 14) 


2. 3 Results 

Calculations were performed for pressure levels of 1, 3, 10 and 30 
atm, for powers of 10 kW and 5 MW. The temperature T q of the 
incoming flow was taken as 300K. The incoming enthalpy for 
referenced to molecular hydrogen at OK as the zero of enthalpy, is 

h = 1.4466 x 10 4 T (J/kg) 
o ° 


. 12 . 


The gas constant for molecular hydrogen is = 4133 J/kgK 

The radiation loss term P T requires specification of the constants 
C and C , Eq. (D-l). They are functions of pressure p and plasma 
size R. # However, they have only been worked out for certain sizes and 
pressures, as described in Appendix D. Calculations could be carried out 
for other values of R. and p if required. Therefore, the values used in 

X o 

these calculations were the ones listed in the following table: 


P (kW) 

p (atm) 
0 



10 

l. 3 

2 

2.4 

10 

10,30 

1. 3 

2.4 

5000 

1,3 

2 

0.4 

5000 

10,30 

2. 2 

0.4 


The initial temperature T^ was usually chosen as 10, 000K, though 
some check runs were made at lower values, and for some low 
intensity cases lower values were used, because the peak temperature 
was near to or lower than 10, 000K. 

The integrations were performed with a s’ o size in x which 
permitted the temperature to change by a maximum of 5% in one step, 
but there was also a specified upper limit on step size (DXMAX) to maintain 
accuracy. This varied with pressure, since the thickness of the wave was 
mainly determined by the absorption length, which decreases with increasing 
pressure. 

An important property of the wave is the amount of laser energy 
absorbed. For low laser intensities, the gas is not heated to a high 
enough temperature to absorb all the energy, which defeats the purpose 
of the wave in the propulsion applications of interest here. A useful 
measure of energy absorption is the percentage of I absorbed by the 
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wave when the temperature has reached the minimum value T on 

mm 

the m^ ^ curve shown in Fig. 2. 1. Very little additional absorption 
can occur after this point. 

Table 2. 1 gives the results of the calculations for the one-dimension- 
al LSC wave model, in terms of the input parameters p Q , I q , R., T^, 

DXMAX, and the output parameters T , T . , m , u ; and % 

’ r r max min A e 

I Absorbed. Here u is the incoming gas velocity, which is related to 
o o 

the mass flow per unit area and the incoming gas density by (2. 2) as 

u = m /0 (2. 15) 

o A o 



v . 

f.i 


X . 


U: 

n ■ ■■ 




The first point to note in Table 2. 1 is that the value of T. does not 
affect the results as long as it is well below the peak temperature, and 
below the temperature at which laser absorption and radiation loss 
occur. This can be seen by comparing Cases 5* and 9, which show 
almost the same results for values of T^ = 10,000 K and 5000 K. A 
similar comparison is made for cases 31 and 32 for a 1000 K change in 
T.. This result is to be expected, since the condition (2. 11) applied 

l 

at the initial station is an exact integral of the energy equation (2. 10) if 
the absorption and loss terms are neglected. The only difference that 
could occur for two different choices of T. is a displacement of the 
origin of x , which is located where T = T^. However, if the origin is 
shifted to the same value of T for two cases with different values of 
T., the profiles cannot be distinguished. 

An example of a temperature profile is shown in Fig. 2. 2 from 
Case 3. The wave at this low pressure is more than 11 cm thick from 
the 10, 000 K point. The divergence between the upper and lower limit 
mass flows begins at 10. 6 cm. The upper limit profile begins to rise 
at 11.6 cm, 15,000 K, while the lower limit profile drops rapidly after 
14,600 K at 11.6 cm, crossing T = 0 before 12 cm. The mass flows 
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exceeded 30, 000K, which is outside the range of conductivities available in the 
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for these u PP er and lower limit profiles differ by about l. 44E-4 
percent. The temperature rise at the front of the wave is quite steep, 
rising from 10, 000 K to 15, 000 K in about 0. 1 cm. This rise is caused 
by heating of the oncoming gas through the mechanisms of thermal and 
radiative energy transport. Only about 4% of the laser energy has been 
absorbed by the 15,000 K point. The major absorption occurs between 
there and the maximum temperature point, 17, 500 K at 2. 6 cm. By 
this point 73% of the laser energy has been absorbed. The high temp- 
erature region is quite broad, with more than 4. 5 cm above 17,000 K. 
After that, there is a gradual fall as the radiation losses exceed the 
heating due to the combined effects of further absorption and energy 
transport by conduction and convection. The absorption is 99% complete 
at 8. 6 cm, where the temperature has fallen to 16, 000 K. 

The very steep rise in temperature at the leading edge of the 
wave predicted by the present model should be viewed with some 
caution. The radiation conduction approximation is not really valid 
in regions of such steep temperature changes, so the details of the 
wave there may not be described correctly. However, energy is 
properly conserved in the model, so the details of the leading edge 
profile should not affect the subsequent features of the wave. In fact, 
observations of the temperature profile of LSC waves in air do show 
temperature gradients of the order of that shown in Fig. 2. 2 (Refs. 

8 aod 9) so the present model is at least qualitatively in agreement with 
observation, even if the exact profile is in some doubt. 


■■2 0 - 
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The general features of the temperature profile shown in Fig. 2. 2 
are typical of all the cases calculated. They all rise sharply, have 
a peak which is broad relative to their total thickness, and then fall 
gradually. The most notable differences are related to pressure and 
intensity. The higher the intensity, the higher the maximum tempera- 
ture, almost independent of pressure. The pressure, on the other 
hand, strongly effects the thickness of the wave, because the absorption 
coefficient controls the rate of absorption of laser energy, and this 
coefficient is strongly pressure-dependent, see Appendix B. The higher 
the pressure, the thinner the wave. For example. Case 23* is the same 
as Case 3 but at 30 atm The peak temperature, 16,900 K, is reached 
0. 032 cm past the 10,000 K point, instead of 2.6 cm past, as with Case 
3. The 15,000 K point is at 0. 092 cm rather than 11.6 cm. The factor 
of 30 in pressure compresses the thickness scale of the wave by a 

factor of 100. 

Perhaps the most important result of LSC wave calculations is the 
relation between intensity and mass flow. For a given value of I q , the 
resulting value of m calculated provides the wave speed U q at which 
the incoming gas must flow, according to (2. 15). For the P = 5 MW 
case a plot of I vs rh is shown in Fig. 2. 3a, while a similar plot 
for the P = 10 kW case is shown in Fig. 2. 3b. On these figures there 
is a line for each pressure, 1, 3, 10, 30 atm. The use of m A rather 
than u permits the curves for all pressure levels to appear on the 
same graph. To give a feeling for the velocities u q , the top horizontal 
scale is the velocity for 3 atm. Other velocities can be obtained easily, 
since they are inversely proportional to the pressure. The velocities 
are listed in Table 2. 1. Also marked on the figure is the percent of 
laser energy absorbed in the wave, if it is less than 100%. 
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From Fig. 2. 3we see that the pressure effect on m is not large 
at any 1^, varying from a factor of 2 for 5 MW at 10 W/cm to a factor 

of 25 for 10 kW at 3E4 W/cm , when the pressure changes by a factor of 
30. 

At the higher intensities, m is approximately proportional to I , on the 
log- log scale, and this dependence persists, for 5 MW, down to 3E4 W/cm , 
which is the intensity region in which more than 90% of the laser energy is 
absorbed. Below 3E4 W/cm the mass flow changes more rapidly 
than the first power of intensity. For 10 kW, the proportionality 
also holds where more than 90% of the laser energy is absorbed. The 
curves bend toward lower m at higher intensities for the 10 kW case 

A 

than for the 5 MW case. 

If attempts are made to calculate m for lower and lower intensities. 

fk 9 

eventually the values of 8®^ very low and finally, no positive value 
will provide a solution which gives an upper limit type of curve for the 

temperature. This is indicative of the threshold phenomenon which has 
been found previously for LSC waves in air. There is some mass flow 
below which the wave cannot be sustained. No attempt was made to 
find this threshold in the present 1-D model for several reasons. First, 
losses may become important near the threshold* and the losses due to 
radial heat conduction are not included here* They are included 
approximately in the model described in section 4. 1 and are there shown 
to have a large effect on the threshold in some cases. Second, the mass 

flows near threshold are very low, and it is felt that values below 10 " 4 

. 2 

g/cm s are not very interesting. Third, the fraction of laser energy 

absorbed near threshold may be very small. Since the main thrust of 

the present work is toward heating the gas, there is no real interest in 

dealing with cases where the laser energy mostly penetrates the wave 

rather than being absorbed. Fourth, at low intensities the T in the wave 

max 
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are near 10, OOOK, and at these temperatures the radiation conduction 
model is not a good approximation to the radiation energy transport, so 
the present model is not a satisfactory representation of the physical 
situation. 

The l-D model of an LSC wave described here can predict the wave 
structure and properties. The most uncertain aspect of this l-D model 
is the choice of constants C 3 andC 2 in the radiation loss expression. 
They should be functions of both pressure and radius R. , but such an 

i 

extensive program of calculation to obtain them was beyond the scope of 
the present work. It is expected that more accurate values of these 
constants will produce only small quantitative changes in the results 
presented here. More serious limitations to the present results are 
imposed by the l-D nature of the model. The effect of several 2-D 
phenomena are discussed in Section 4. For sufficiently large size 

plasmas, this l-D model is probably applicable. 

Whether this model realistically describes LSC waves in hydrogen 
can only be determined by comparison with experimental measurements 
of such features as and m^. There are no experimental results 

for LSC waves in hydrogen available at the present time. 



3. COUPLED LSC WAVE. STREAMTUBE FLOW 


3. 1 Basic Equations 

The flow of the gas in the streamtube brings in two new phenomena 
not included in the 1-D LSC wave model. These are area changes and 
pressure changes. The purpose of the area change is to accelerate the 
gas, in order recover maximum momentum for purposes of propulsion. 

With these new phenomena, all the equations of motion change. The 
mass conservation equation now includes the cross-sectional area A: 


. • 2 
puA = m = constant, A = ttR (3. 1) 

The momentum conservation equation expresses the balance of pressure 
and inertia forces: 


dp du 

+ P U S " 0 


(3.2) 


The energy conservation equation must now include the area terms, so 
(2. 1) is modified to 


5 <” H > + I < AI > + Ap x 


S [ A <\= +»*>£] 


(3.3) 


Here we have again used the conduction approximation for the radiation 
energy transport term. 

The inclusion of the terms on the right in a streamtube calculation 
deserves some comment. Ordinarily, streamtube calculations do not 
include these axial conduction energy transport terms, and in fact they were 


- 27 - 


ORIGNM-t -AjGBO 

OF POOR Q UAU ™ 




not included in the energy equation of Ref. 1. However, for LSC waves 
the temperatures are sufficiently high so that the radiation conduction coef- 
ficient is large, making the axial radiation transport significant even 
R 

though the temperature gradients are not large. Once tae radiation transport 
is included, the thermal conduction term can also be carried along with very 
little additional effort, and it also may be significant at the higher temperatures. 
The inclusion of these transport terms makes the energy equation second 
order in T, rather than first order, as it was in Ref. 1, so an additional 
boundary condition is needed. 

The reduction of (3. 3) to an equation in T and p is straightforward. 

Since the gas enthalpy is a function of these two variables, its derivative 


dh = c dT - c dp 
P T 


c = (dh/dT) 

P P 


c T = (5h/3p) T 


With m constant, the first term of (3. 3) can be written 


. / dh du \ 

m \ dx + U dx / 


• / dT dp du 

m [c — + c_ ” + u — 
1 p dx T dx dx 


(3.4) 


The laser absorption term again follows from the energy transport equation 


dAI/dx = k^AI 


AI = Vo*'' ’ T = / k L dX 




where A is the cross-sectional area of the LSC wave plasma. Then 




d Al/dx = k. A. I e' 
L 1 o 


We also note that this nozzle flow starts at the back of the LSC wave so that 


m = m A A i* Then we can U8e < 3 * 4 ) and (3. 6) in (3. 3) and divide by the 
constant m to get 


d£ du d [/ X C + X R \dT k L J o C T 

°T dx U dx " pu dx || p u ) dxj + m A 


Since p can be eliminated by the equation of state (2. 7), we have three 
equations (3. 1), (3. 2) and (3. 6) for the four variables p, T, u and A. When 
one of them is specified as a function of x, the other three can be determined. 
The conventional way to solve this problem is to specify a streamtube shape 
A(x) and solve for p, T, and u. However, as pointed out in Ref. 1, this 
method leads to a singular point when u equals the local sound speed, which 
makes the integration a difficult numerical problem. A much simpler inte- 
gration can be done if u(x) is specified and the equations are solved for p, T 
and A (or R = (A/tt)*^). This was the method used in Ref. 1 and it will also 


be used here. 


To this end, (3. 2) is used to write the left side of (3. 6) as 


dT dp du dT 

c ~ + c_ -r 6 + u — = c T" + (1 
p dx T dx dx p dx 


\ du 
PC T> U S 


The coefficient of the second term is given in Eq. (A-l) of Appendix A as 


1 - p c T = - (b In p In T) 


and differentiation of the equation of state (2. 7) gives 
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1 - p c„ = 1 + (3 In Z/3 In T) 

1 P 

Then the energy equation (3. 6) is 



(3.7) 


Equations (3. 1), (3. 2) and (3. 7) are now three rather simple equations 
for p, T and R when u(x) is given. To them must be appended the definition 
of t, (3. 5). For computation they are written as a first order system by 
defining 


v - X C * X R dT 
pu dx 


(3.8) 


so the system becomes 


dT _ pu V 
dx + 


(3. 9a) 


dV 

dx 


c o u V 




+ X 


1 + 


R 


3 In Z 
3 In T 


1 du + P T 

u — + - 

J d* P u 


k. I 
L o 

m. 


-T 


(3.9b) 


dp/dx = - pu du/dx (3.9c) 

dt/dx = k. (3. 9d) 

ia 

R/R. = (m A /pu) 1/2 (3.9e) 
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The density p is given by the state equation (2. 7) and the thermodynamic 
derivative in (3,9b) is given in Appendix A, Eqs. (A-19a) and (A-19b). 

The initial conditions on most of the variables are straightforward. 

One begins the integration at some value of x, where the values of T, u, 

р, t are known from the LSC wave solution. If we denote the point where 
the streamtube solution is joined to the LSC wave solution by the subscript 

с, then we have 

x = x: T = T , p = p, T = T , u = u 

c C O c c 

One further initial condition, on dT/dx, is needed. The specification of 
that condition is not so straightforward, and it is discussed next in connection 
with the description of the numerical method of solution. 

Notice that we perform the integration in the physical distance variable 
x. In Ref. 1 we used the optical depth variable t instead. This variable has 
some utility when the radiation transport terms are not included, because 
then k does not appear in the energy and momentum equations. It only 

Ij 

serves to define the physical distance scale through (3. 9d). However, in 
the present case, where P T and the transport terms are included, there 
is no such simplification gained by using t, and we chose to use x as the 
independent variable. 

3. 2 Numerical Method 

The joining of the streamtube equations (3. 9) to an LSC wave solution 
requires one more boundary, on dT/dx, in addition to those given in (3. 10). 

If one merely picked a station x c and used the value of dT/dx there, the 
streamtube solution would either rise or fall, just as the LSC wave solution 
did. One needs a degree of freedom for the streamtube integration cor- 
responding to m in the LSC wave integration. This freedom is provided 
by dT/dx. Near the rear of the LSC wave solution the values of dT/dx (or 

-31- 

ORIGINAL PAGE IS 
OF POOR QUALITY 


W, see Eq. (2. 12a) ) are different for the upper and lower limit solutions, 
even though they correspond to solutions with nearly the same value of 
m . Thus, dT/dx is only determined within some range. For example, 
in the case to be discussed here. Case 13', at x = 2. 97 cm, the values of W 
differ by 17%. Thus, if we start a streamtube solution there, one may use 
any value of W between the two limits to obtain the initial value of dT/dx 
(or V, see Eq. (3. 8) ), and still be consistent with the LSC wave solution. 

With this degree of freedom available to us, we may set up a scheme for 
joining the streamtube calculation to an existing LSC wave solution. 

First, we wish to choose a matching point x . It seems desirable to 
choose this as far toward the back of the LSC wave as possible, so that most 
of the gas heating and laser absorption occurs in the constant area wave, 
which acts as a "combustion chamber". We therefore choose the upper 
limit LSC wave solution, and choose a number of points near the end of 
the wave, before the minimum temperature point. At each point we inte- 
grate the streamtube equations, using the value ofV = W/pu = W/rii^ 
appropriate to that point. For smaller values of x, the resulting tempera- 
ture profile will eventually rise. We find the smallest value of x for which 
it rises, and choose this as our matching point x^ 

We then integrate the streamtube equations again from this x.^, but 
this time we allow the initial value of W to vary, and we iterate on it to find 
upper and lower profiles of T, just as we did for m A in the LSC wave solution. 
This integration will progress a certain distance in x, until the T profiles 
diverge, one to fall and one to rise. We then pick a point where T differs 
by only a few degrees on the upper and lower limit curves, start integrating 
there, and again iterate on V to find another pair of limit curves. We con- 
tinue in this fashion, moving toward large x step-by-step for as long as we 
choose, to calculate the conditions in a streamtube matched to any given LSC 
wave solution. 




The results of one sample calculation using this method are given 
below, and they show this to be a feasible scheme of calculation, which 
produces a consistent solution. There may, of course, be other schemes which 
also work. No attempt was made to find the most efficient scheme. 

Before the results are given, we will describe the velocity distribution 
function chosen for our sample case. This function is a required input to 
the streamtube integration. 

3.3 Velocity Distribution 

The form of the function u(x) chosen will determine the shape of the 
streamtube R(x). Several conditions on u(x) are clearly desirable. At x 
it should join smoothly on to the velocity in the LSC wave, where value u C 
and slope (du/dx) c are known. It should increase to a high value since we 
wish to accelerate the flow in the streamtube, and in fact should reach the 
local speed of sound at about the location where the throat is desired. 

There are clearly many functions u(x) which have these properties. 
Whether any given function is satisfactory can only be determined after the 
calculation is performed and the resulting streamtube shape determined to 
be a reasonable one. 

For the sample calculation performed here, the velocity distribution 
was taken in the form 



The first term assures that u and du/dx have the values u and (du/dx) 

at x = x c , providing the power E on the second term is greater than unity. 

The second term provides the increase in u as x increases. Since the slope 

of u in the LSC wave is small at the joining point x , the value of u reached 

c 

at x = x^ will be nearly au^. Thus the parameters a and determine the 
velocity reached at some chosen station x^. 
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The slope of this expression is 


(3. 12) 



For the particular case calculated here, the joining point x was about 

c 

3 cm. The point x^ was chosen as 5 cm, and a was taken as 10 so that 
at x = 5 cm, the velocity would be approximately 10 times u . The 
exponent E was taken as 2. Changes in a, E, and x^ clearly will alter the 
length of the streamtube to the throat, as will its size and shape. Experience 
with the calculations should permit choosing the parameters so as to obtain 
any desired streamtube shape of reasonable form. 

3. 4 Results 

The case used to illustrate the coupling calculation was Case 13', 

2 

p =3 atm, I = 3E5 W/cm , F = 5 MW. This was chosen as close to one 
o o 

of the cases calculated in Ref. 1. 

The calculation was completed in 5 steps using the velocity distribution 
function described above. The course of the calculation is given in the following 
Table. The value of x at the beginning of each step is given, starting with 

x . the matching point to the LSC wave. The values of T , T , W and W 

Q U i-/ U 

on the upper and lower curves of the previous step are also given, as well 
as the value of W found by iteration for the step starting at x. 


x(cm) 

T u (K > 

t l (K) 

W y x 10‘ 6 

W T x 10' 6 
L 

W x 10" 6 

2.9683 

17, 129 

17, 098 

-47. 37521 

-56.84703 

-50.48019 

4. 2183 

16, 509 

16 , 508 

-2. 251405 

-2. 457159 

-2. 294881 

4. 6983 

16,333 

16, 328 

-1. 638167 

-3. 242737 

-1. 814397 

5. 0183 

16 , 168 

16 , 1 67 

-1. 864034 

-2. 245204 

-1. 984559 

5. 2883 

15,982 

15, 980 

-1. 817534 

-2. 539966 

-2. 341516 
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It is seen that the step sizes decrease as we approach the streamtube 
throat. The last step passed through the throat before the temperature 
began to rise, and the calculation was terminated with this step. The throat 


conditions were 


x = 5. 5583 cm, T = 15, 740 K, h = 1. 308E9 J/kg, p = 1.77 atm 


u = 15. 07 km/s, R/R^ = 0. 2943, R = 0.677 cm, M = 1.08 


The laser energy was completely absorbed, and in fact, at x^ = 2. 9683 cm 
only 4. 9E-3 percent of the laser was not absorbed, so the streamtube calcu- 
lation could have proceeded with the laser absorption term. 

Notice that the value of W resulting from the iteration at each step 
was between the upper and lower values for that location in the previous 
step, which testifies to the consistency of the calculation scheme. The 
uncertainty in W can be used to make the calculation progress downstream 


in a consistent manner. 


Note that the Mach number M at the physical throat need not be unity 
when heat addition is occurring. The pressure is down from 3 to 1.77 atm, 
which is a little higher than the perfect gas (y = 1.4) value of 1.58 atm, be- 
cause of real gas effects. In connection with the pressure it should be noted 


that X„ was only tabulated at specific pressures, and the present calculations 


were performed by using the value at the nearest pressure, rather than by 
interpolating. It would probably be somewhat more accurate to interpolate 
in pressure, at the expense of some complication in the program. 

Plots of T (x), u(x) and R(x)/R i are given in Fig. 3. 1. The tempera- 
ture profile shows the characteristic rapid rise, peak and slow decay. 

The maximum is at 19, 000 K at 0. 5 cm. The tic marks on this curve, 
and the others, indicate the starting points for the integrations, as listed 
in the above table. In the streamtube portion of the flow, the temperature 
only drops from 17, 129 K to 15, 740 K, 
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5 MW, 











The velocity profile shows a rapid rise to the station of maximum 
temperature, and is then nearly flat, with a slow decrease, during the LSC 
wave. This corresponds to the density showing a slow increase, since 
p u is constant. This is consistent with the temperature drop of 1400 K 
in this region, since p is constant, and the gas is fully ionized, so pT 
is nearly constant. During the streamtube portion of the flow the velocity 
increases rapidly according to the relation given in (3. 11) and the constants 
specified below (3. 12). It reaches 10 u at x = 5 cm, as specified and 
continues to increase to the throat. 

The radius ratio remains at unity through the LSC wave, and then 
decreases in the streamtube region, reaching the throat value of 0. 2943 at 
x « 5.56 cm. This shape could have been altered by changing the velocity 
distribution parameters. If x f had been made larger, the streamtube would 
have been longer, probably with nearly the same throat size, but with a lower 
temperature due to the additional radiation losses brought about by the 
increased length. 

For this calculation one can easily estimate the total radiation losses 
by finding the energy contained in the flowing gas at the throat. With the 

2 

values of h and u given above, the total gas enthalpy is 1 308E9 + 0. 5 (1. 507E4) = 

2 

1.42E9 J/kg. The mass flow per unit area in the wave is 1. 228 kg/m s and 
the radius of the plasma in the wave is 0. 023 m, so the total mass flow is 
1. 228 tt ( 0 . 023) 2 = 2. 04E-3 kg/s. The total power flow through the throat is thus 
2.90 MW, which is 58% of the input power. Thus 42% of the power has been 


lost in radiation. 

This calculation may be compared with one of the cases in Ref. 1. 
There, the presence of the LSC wave was not considered, except to estimate 


the mass flow to support the wave by a very crude calculation of the wave 

2 

structure. The result, for P = 5 MW, p = 3 atm, I = 3. 67E5 W/cm , was 

2 o o 

m. = 0.3316 g/cm a a value plotted as the circle on Fig. 2.3a marked Ref. 1, 
A 2 

It is a factor of 2. 7 higher than the value of 0. 1228 g/cm s calculated in the 
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present work for I q = 3E5 W/cm . With this high value used, the calculation 

of the streamtube in Ref. 1 was carried out by using the equations of this 
section, but without the energy transport terms on the right of (3. 3). 

The calculation was started at 12, 000 K, with the laser energy reduced 
by the amount of energy necessary to heat the assumed mass flow to 
this temperature. The resulting temperature distribution is shown 
in Fig. 3. 2, taken from Fig. 7. 2a of Ref. 1. Without the presence of 
the LSC wave, and with nearly 3 times the mass flow, this profile 
shows a much different character than that of the present work. Fig. 

3. 1. It rises gradually to a peak of only 14,900K, and drops to 14, 300K 
at the throat. The lower temperature is largely a result of the higher 
mass flow at nearly the same laser intensity. But lower temperature 
reduces the radiative losses, as does the shorter length of the flow, 3. 2 cm 
instead of 5, 5. The result is that the calculation in Ref. 1 of the radiation 

losses showed them to be only about 5% of the input power, rather than 
the 42% in the present more accurate calculation. 

This comparison emphasizes the necessity of having a correct 
calculation of the mass flux through the LSC wave in order to properly 
obtain the temperature levels in the flow, and the resulting losses. 

Based on this example, the radiation loss calculations in Ref. 1 must 
be considered as doubtful. In order to assess them more accurately, 
streamtube calculations by the present method should be coupled to 
the LSC wave calculations described in Section 2. The single sample 
calculation described here shows that the present method can be used 
for this coupling. 

The large losses found in the present calculation are attributable to 

the high temperature level in the gas. This in turn is caused by the necessity 
of having a significant degree of ionization, so that the absorption mechanism 
of inverse bremsstrahlung can operate. Ionization in hydrogen requires 
temperatures in excess of 10, 000K. It would be advantageous to operate 
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at lower temperatures to reduce the losses. This could be achieved by using 
an easily ionizable seed in the hydrogen, to produce electrons at lower tem<- 
perature. However, the type and amount of seed must be chosen with care. 
The weight of the seed will reduce the specific impulse of the gas, and so 
will the reduced temperature. It would be desirable to conduct a study 
similar to the present one, but using seeded hydrogen, to determine the 
effect on mass flux, specific impulse, and losses, of various types and 


amounts of seed. 
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4. TWO DIMENSIONAL EFFECTS 


J - V- 


4. 1 Heat Losses 

One of the two-dimensional effects not included fully in the one- 
dimensional model is transverse (radial) heat loss from the plasma cylinder. 
The only loss of this type included is the radiation loss from small cylinders, 
where the plasma is transparent. For larger cylinders, the radial radiation 
loss also includes black radiation from the edges of the cylinder due to the 
radial temperature gradient, which is not accounted for in the one- dimensional 
model, as explained in Appendix D. And for all cylinder sizes, the transverse 
conductive loss has not been included. In this subsection we describe an 
estimate of these losses and show their effect on the m A vs I q relation. 

To account for radial conduction losses a term 



has to be added to the right side of Eq. (2. 1), where r is the cylindrical 
radial coordinate. To include this in the framework of the one -dimensional 
model one needs to express it as a function of x only. Let us assume q yc 
depends on x only. We may then integrate (4. 1) twice on r between the center 
of the cylinder r = 0 and the edge, R., to get 


q 


rc 



x c dT 


(4. 2) 


where T is the temperature at the edge of the cylinder, 
o 
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T 


# 


An expression of the form (4. 2) provides an estimate of the heat loss 
term which is to be added to the right side of (2. 1) to account for transverse 
conduction loss. The value to use for T is not certain since we do not 
know the radial temperature profile. Since most of the integral is at high 
temperature, the value of the low edge temperature will not have much effect 
on q rc , so we take T g = 300 K, which gives the integral its largest value. 

To compensate for this we could reduce the factor 4 somewhat. In fact, 

(4. 2) has already been suggested by Raizer in Ref. 2 to estimate the conduction 
losses, thougn ne used 2. 9 instead of 4. We adopted the value 3 in our 
calculations. 

In addition to the conductive loss, there is also the radiation loss for 
large plasma cylinders, as mentioned above. If we use the radiation conduction 
approximation for this term, it is also of the form of (4. 2), with \ R instead 
of \ c . 

To use these estimates for assessing the effect of radiation losses 
we simply add to the right side of (2. 12b) the terms 


X c dT + 


"lr f 

R 2 J 

i 300 


X R d T 


(4.3) 


The first integral is tabulated in Table C2, Appendix C as 1^ . It is a 
function of T and quadratic interpolation was used between the tabulated values 
of T. The second integral is tabulated in Table Dl of Appendix D as L and 

*R 

was used in the same way. 

Calculations were made including these loss terms by the same scheme 
described in Section 2, 2, using the values 


'lc 

= 3, 

°LR * 3 

for 

P = 5 MW 

LC 

= 3, 

o 

II 

U 

for 

P = 10 kW 
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The results are shown in Fig. 4. 1, which is a plot of I vs m A . The lines 

are the results for no sideways losses from Fig. 2,3, while the symbols 

show the results including losses. For the 5 MW power. Fig. 4. la, there 

is very little effect of these losses, even at the highest intensities, where 

R. is the smallest and the temperature the highest, which make the loss 

terms the largest. But for the 10 kW power. Fig. 4. lb, these losses show 

a large effect because of the small sizes involved. For 3 atm, solutions could 

be obtained only at I = 1E6 and 3E5 W/cm , but there was no solution which 

°2 

has ™A > lE " 4 S/c m 8 at 1 = 2E5. That intensity is marked as the 3 atm 
threshold on Fig. 4. lb. Similar thresholds are marked for 10 and 30 atm. 

For P Q = 1 atm » n ° »^ch solution was found even at 1E6 W/cm 2 . Thus, 
when the losses are included, thresholds are found much above the intensities 
at which they would have been found without losses, for the 10 kW power. 

Where solutions were found, they occur at lower mass flow for a given 
intensity. The reductions in mass flow are quite small at the higher pressure, 
but are large at the lower pressure; at 3 atm the reduction is about a factor 
of 2. This large effect at low pressure is caused by the large length of the 
wave, which allows more area to contribute to the radial heat loss. 

This limited set of calculations, which is tabulated in Table 4. 1, shows 
the importance of including the transverse loss when considering small sizes 
and low pressures. The effects found are believed to be qualitatively accurate, 
but probably not quantitatively precise. The expressions (4. 3) are of the 
correct form, and have the correct dependence on radius, but they are one- 
dimensional approximations to a two-dimensional phenomenon, and contain 
coefficients which can only be estimated roughly. A more precise eval- 
uation of the transverse heat loss due to thermal condition and black radiation 
requires a two-dimensional model of the wave. 

i 

& 
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Fig. 4. la Intensity vs Mass Flow for a Steady I. SC Wave with Radial Toss — P = 5 MY. 







TABLE 4. 1 

RESULTS OF LSC WAVE MODEL WITH RADIAL HEAT LOSS 
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4. 2 Beam Convergence 

The convergence of the laser beam provides another two-dimensional 
effect which can change the relation between intensity and mass flux in an 


LSC wave. If the cross-sectional area of the beam A^, varies with x, then 


the intensity I will vary not only because of absorption, but also because of 
the area change. The conservation of total power in the laser beam requires 


that IA change by absorption, according to 

IJ 


d IA 


B 


dx 




(4.4) 


this is the same relation as (3. 5), except here the appropriate area is the 

beam area A ; 

U 

before, namely. 


beam area A instead of the flow area A, The solution is the same as 

Jj 


where the subscript i refers to the location x = 0, 

In the one-dimensional model of the wave, we are concerned with con- 
ditions near the centerline of the beam, and the laser absorption term in the 
energy equation is the same as in (2. 1) and (2.4), namely. 



LA. = I. A e 
B i Bi 


T = 


7 


k dx 

.Lf 


(4.5) 


“ ■ = - k I 
dx L 


(4.6) 


Use of (4. 5) then shows this term to be 


dl 

dx 


A 

i -T t Bi 
k T e I. -r — 
L i A fi 
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fl 


(4.7) 


This shows that the modification required for consideration of beam conver- 
gence is to insert the area ratio factor Ag./Ag i„ the laser absorption term 
This requires some specification oi the beam shape. 

A suitable expression for this shape is found for example, in Ref. 10, 

where the intensity is given as a function of radius r and distance from the ’ 
focal spot z, as 



The shape of the beam is Gaussian in the radial direction, and parabolic along 
the beam axis, as defined by the function cju(z). The minimum area is at 
There is a size parameter, represented by or m , which also 

depend, on the laser wave length A .ketch of the .hape°i. ahown in 
Fig. 4. 2. 

The parameter is determined by the geometry of the particular lens 
system used in any given case. Suppose the lens is a distance z = -L, from 
the focal plane, as shown in Fig. 4.2. The corresponding value of uj is de- 
noted as At this value of z, (4. 8) becomes a quadratic equation for 

oi whose solution is 
o 
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Fig. 4. 2 Schematic of Converging Laser Beam 
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Since the wavelex^fth of interest here is small, 10. 6 nm, the radical may be 
expanded to give 





(4.9) 




i 


m 


■ 




the quantity L/u) T is the ratio of the distance from the lens to the focal plane 

J-i 

to the size of the beam on the lens, which is the f- number of the lens, Thus, 
the parameter z is related to the wavelength and a known property of the 
lens system. 

The beam area can be defined as extending to some fixed point on the 

2 

Gaussian decay curve, that is to r/u) = R_/uj = constant. Thus R is pre- 
2 B B 

portional to id . Suppose the beginning of the wave, denoted by subscript i, 
is at a distance Xp from the focal plane as shown in Fig. 4. 2. Then R Bi is 
at z = - x p> and any other location in the wave is at z = x - Xp. The area 
ratio is then 



2 / x 2 

z + (x-x ) 

o J 


2 2 
z o + *F 


(4. 10) 


This defines the new term in the energy equation in terms of the two physical 
parameters W^ L and Xp. The former is the lens f-number as mentioned 
above, and the latter is the distance from the focal plane to the wave 
origin. Notice that the intensity I. in Eq. (4. 7). is the value at x = 0, 
where the wave originates. 

With the factor A ./A inserted in the absorption term on the right 
Bi B 

of (2. 12b), and (4. 10) used for this ratio, the LSC wave energy equation 
may be solved as before, once L/uj^ and Xp are specified. The case of a 
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* 


constant beam area, considered in Section 2, corresponds to L/cu very 

L 

large compared to x F , since then (4. 10) shows the area ratio to be unity. 

As an example of the effect of beam convergence, the calculation 

of Case 103 was re-done. This is the case of P = 10 kW, p =3 atm 
2 o ’ 

I. = 3E5 W/cm . The value of L/ w ^ was taken as 30, which is a high 

f- number. The value of x p was estimated by taking the edge of the beam 
at the e' level, so that 2 r|./ iju 2 (x f ) = 2. With R fi . determined from the 
power and intensity by P = ttR 2 .I. ( with 'ju(x F ) related to z by (4. 8), and 
with z q from (4. 9), the value of x is 3. 075 cm. The LSC wave calcula- 
tion then gives m^ = 0. 0845 g/cm s, u^ = 344 cm/s, compared to 0. 0730 
and 297 for a parallel beam (L/u, L «), Case 103. Thus this rather mildly 
convergent beam, which has a focal spot radius of 0. 01 cm, increases the 
mass flux by 16%. The location at which 99% of the beam is absorbed is 
1 . 44 cm, well before the focal spot location of x_ = 3. 075 cm. The maxi- 
mum temperature reached is 20, 100K, an increase from 19, 000K for the 
parallel beam. 

A convergent beam behaves as a higher intensity beam, relative to a 
parallel beam of the same initial intensity, and necessitates a higher inflow 
velocity. No other calculations of waves in a converging beam were made. 
However, they can easily be performed when the geometric parameters 
L/o)^ and x F are specified. 

A converging beam has the desirable feature of providing static sta- 
bility for the position of the LSC wave, because of the increase in mass flux 
with intensity characteristic of LSC waves. If a wave in a convergent beam 
is displaced upstream to a lower intensity position, the mass flux to keep it 
Steady there will be lower than the mass flux actually occurring, and it will 
be blown back downstream. If it is displaced downstream to a position of higher 
intensity, the mass flux required at that position will be higher than that occur- 
ring, and it will try to increase the mass flux by moving upstream. In the 
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divergent portion of the beam, downstream of the focal spot, the same argu- 
ment shows the wave position to be statically unstable. The question of 
dynamic stability needs further investigation. 

4. 3 Transverse Velocit y 

In the LSC wave experiments done in air, it has been recognized that 
transverse velocities may be occurring. These experiments are not confined 
in the radial dimension, so the beam only fills a small part of gas volume. It 
is then possible for some of the oncoming flow to turn radially and flow around 
the wave rather than through it. If this actually occurred, it would have the 
effect of making the experimentally measured wave speeds higher than the ones 
calculated by a one-dimensional model. This can be seen as follows. (Fig. 4. 3) 
If some of the flow goes around the wave then the streamtube of flow going through 
the wave has a larger area when it enters the wave than it does far upstream. 

Thus its axial speed is lower as it enters the wave than it is far upstream. 

In the one-dimensional model, we calculate the speed as it enters the wave. 

But the experiments measure the speed far upstream (wave speed relative to 
the undisturbed flow). Thus the experimental speeds could be higher than the 
speed of the gas entering the wave. 

The measurements made on LSC wave speeds in air do, in fact, show 
them to be higher than the one -dimensional predictions by about a factor 
of 3, Ref. 4. They also show evidence of two-dimensional flow, which is 
consistent with the above explanation. Therefore, in unconfined flows, the 
one -dimensional predictions of wave speed could be too low because of the 
neglect of transverse flow effects. 

The only attempt made to calculate the effect of transverse velocity- 
on LSC wave structure that is known to the authors is that of Ref. 11, 
which was quite crude. They estimated that the radial flow might cause 
the inflow far from the wave to be 1. 5 to 2. 25 times the value predicted 
by two-dimensional theory. 


- 53 - 







In th. propulsion application of inter.at in the pre.ant work, th. flow 
occurs in a confined channel, with the laeer beam enviaioned a. heating a 
major fraction of the gas flow, and with walls bounding the flow in the 
radial direction. In this case the flow does not have as much freedom to 
develop radial motion as in the unconfined case. It would appear that the 
effects of transverse velocity should be much less than in the unconfined case. 
Unfortunately, no experiment, on LSC wave, in confined channels have been 
reported, so empirical evidence for this view is not as yet available. 

Considering the difficulty of constructing a model including transverse 
velocity, and lack of evidence that it is an important factor in confined LSC 
waves, it was not considered useful within the limited scope of the present 
work to pursue the subject further. 
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5. CONCLUSIONS 


Y/e have constructed a one-dimensional model of a laser supported 
combustion (LSC) wave in hydrogen. The pressure in the wave, and the 
mass flux per unit area, are constants. The temperature profile is deter- 
mined by the energy conservation equation as a balance between convection, 
absorption of laser energy, radiation loss to the surroundings, thermal con- 
duction, and radiation transport of energy from the hot part of the wave to 
the cool part. It has been found that for the conditions of interest the radia- 
tion transport may be treated in the radiation conduction approximation. 

To obtain the wave structure for a given pressure and laser intensity, 
a numerical iteration must be performed to find the mass flux. Calculations 
have been made for pressures of 1,3, 10 and 30 atm, for 10. 6 urn laser 
intensities from 10 4 to 10 6 W/cm 2 . The resulting temperature profiles all 
show a rapid rise due to heating by conduction and radiation, a high peak 
due to absorption of laser energy, and a slow decay dominated by radiation 
losses. The peak temperatures are high because the laser absorption can 
only occur by inverse bremsstrahlung, and this mechanism requires elec- 
trons, which in turn requires high temperatures in hydrogen. At the highest 
intensities and pressures, the maximum temperature ranges between 
20, 000 to 30, 000K. 

The mass flux at each pressure and intensity determines the velocity 
with which the gas must flow into the wave to keep it stationary. The mass 
flux increases with pressure at a given intensity. At low intensities the 
mass flux becomes very small, dropping off toward zero at some threshold 
value of intensity. This threshold value is strongly affected by transverse 
losses. 
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In a laser- heated hydrogen rocket, it is necessary to expand the gas 
issuing from the LSC wave through a converging-diverging nozzle. The nozzle 
was modeled by a variable -area, variable -pressure streamtube flow, by 
adding appropriate terms to the LSC wave energy, mass and momentum con- 
servation equations. A computational scheme was developed for the stream- 
tube flow, using the conditions near the back of the LSC wave as initial con- 
ditions. One case was calculated to prove the usefulness of the scheme, 
and the results compared with a calculation from Ref. 1. In that study, the 
wave structure was not calculated, so the mass flux to hold the wave station- 
ary was only estimated. The present result showed that the mass flux found 
by solving the wave structure was nearly a factor of 3 less than that estimated 
in Ret. 1. This causes the maximum temperature to be much higher than 
found m Ref. 1(19, 000K rather than 15, 000K), which in turn makes the radia- 
tion losses much larger. The energy remaining in the gas at the streamtube 
throat is only 58% of the laser energy, rather than the 95% found in Ref. 1. 

The present more complete calculations emphasize the high tempera- 
tures and high attendant losses which come from the combination of 10. 6 (jm 
laser radiation and hydrogen working fluid. They are caused by reliance on 
the electron-dependent mechanism of inverse bremsstrahlung for laser 
energy absorption in a gas where electrons are produced in quantity only 
above 10, 000K. It appears that use of an easily ionizable seedant would be 
advantageous to permit absorption at lower temperatures and reduce losses. 

Several two-dimensional effects were examined. The first was radial 
heat losses, caused by conductive heat transfer to the surroundings and by 
black radiation from the edges of the hot gas. Both were modeled approxi- 
mately to fit into the framework of a one -dimensional model. They are 
found to be inversely proportional to the cross-sectional area of the hot gas. 

For the 5 MW power level, these losses have little effect on the mass flux 
at any given intensity, because this area is large enough to make the radial 
losses small. However, for the 10 kW power, the small areas produce large 
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radial losses, which sharply reduce the mass flux from the value it has when 
these losses are ignored. The results show an effective threshold for the 
wave which is above 10 W/cm for 1 atm, and above 10 4 W/cm 2 at 30 atm, 
with intermediate values between. This reinforces the view expressed in 
Ref. 1 that a laser-heated rocket of this type, in a 10 kW size is a loss- 
dominated device, not scalable to higher powers. 

A second two-dimensional effect considered was the convergence of 
the laser beam, causing an increase in intensity. This was modeled by spec- 
ifying a laser beam shape in the one-dimensional absorption term. One calcu- 
lation showed the expected result that the mass flux increased if the beam 
converged from its initial size at the wave front, instead of remaining parallel. 
An argument was presented showing that a converging beam leads to static 
stability of the wave position; dynamic stability was not investigated. 

In view of experimental observations of its occurrence in unconfined 
flows, the question of transverse flow was addressed briefly. It is specu- 
lated that transverse flows cause the large observed deviations of mass flow 
from one-dimensional theoretical values. However, the propulsion applica- 
tions envisioned as the motivation for the present work will involve flows 
which are confined radially. Therefore, transverse flow is expect *d to be 
less severe than in radially unconfined flows, although there are . > experi- 
ments in such flows to verify this assertion. Considering this uncertainty 
of the importance of transverse flow in propulsion applications, and the 
difficulty of estimating such effects theoretically, no further work on t’.is 
subject was undertaken in the present study. 
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APPENDIX A 


THERMODYNAMIC PROPERTIES OF EQUILIBRIUM HYDROGEN 

We need the derivatives of enthalpy with respect to p and T, the 
equation of state which relates p and T, and the composition of the 
mixture of hydrogen molecules, atoms, ions and electrons. The pressure 
derivative of h can be expressed in terms of the equation of state by a 
classical use of reciprocity relations between the derivatives of h and 
entropy, which are in turn obtained from the second law. These lead to 




(el ■ ( 


5 In p 

6 In T 


(A- 1) 


We characterize the hydrogen as a mixture of perfect gases whose 

components have number densities: molecules n w , atoms n , ions n and 

M A I 

electrons n^_,. We do not expect to get a temperature where other species 

are significant, nor are there multiply ionized atoms, so n = n . As 

I E 

reaction coordinates we use the fraction of molecules dissociated 0 and 
the fraction of atoms ionized, a: 


(n A + n^/2 n x 

5 " "M * (n A + n I> /2 ' “ = n A + n I 
In terms of the original number of particles per unit volume. 


(A-2) 


"o - n M + (n A + n I )/2 


(A-3) 


the number density of the species are 


n M = n o (1 ‘ P) ' n A s 2n o 3 (1 " fl) ' n l = n E = 2n o 3a (A_4) 


and the corresponding partial pressure are found from 


p. = n. kT 
1 l 


(A-5) 


where k is the Boltzmann constant. The partial mass densities are, in 


terms of the particle masses m. , 


P M = n M °V P A = n A 'V P I = n l m i' P E = n I m E (A - 6) 


and the total number and mass densities are 


n = n [1 + 3 (i + 2a)], P = n m w 
° o M 


(A-7) 


Whenever convenient we will use the facts that an excellent approxima- 
tion for hydrogen is m. = m T = m,,/2, 

AIM 

The equation of state is the sum of the partial pressures. 


p = ° 2 R M T » Z = 1 + Ml + 2a), R m = k/m M 


(5 In p/d In T) = - 1 - (d In Z/d In T) 

P p 


(A-8) 


where Z is the compressibility factor and R. , the gas constant for the 
molecules. We have now expressed p in terms of p, T through a and 0 

by (A-8). Since p and T are the dependent variables, we need a, 8 in 
terms of T, p. 
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Since we deal with hydrogen in thermochemical equilibrium, a and 
0 are determined by the Law of Mass Action in terms of the thermodynamic 
variables of the mixture. There are two reactions occuring, dissociation 
and ionization: 

H 2 # 2H H £ H + + E~ 

For each the equilibrium constant provides a relation between the 
partial pressures and the partition functions of the species, involving the 
heat of reaction. These relations are 


M 


TTm 


3/2 


(kT) 5/2 Q 


elA 


- h D m M /kT 


Q Q Q Q ... 

c v r elM 


(A-9a) 


P I P E 


2rrm. 


3/2 


(kT) 


5/2 Q ell °elE 


Q 


-h m /kT 
e (A-9b) 


elA 


where h„ is Planck's constant, Q and Q are the vibrational and rotational 
P v r 

partition functions of H , Q is a partition function reflecting vibration- 

2 C 

rotation coupling, Q^. are the electronic partition functions of the species, 

h ° is the heat of dissociation of a molecule per unit mass and h ° is 
D ■*••*-* 

the heat of ionization of an atom per unit mass. From (A-4), (A-5), (A-7) 
and (A-8) the partial pressure/ratios are 


M 


, kT . ^ 2 a -.) 2 i 

1-0 o 1-0 Z 


(A- 10a) 


. aiaL. n kT 

P A 20 ( l - a) o 
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The electronic partition function for the molecule may be replaced by its 

ground- state statistical weight 1, since its first excited state lies very 

high (132* 000 K). For the electron Q^ E = 2 because of the two electron 

spin states, and for the ion, there are no internal states, Q =1. 

ell 

Using (A- 10) in (A-9), with Z from (A-8), yields 


- P 2 (l-tt ) 2 


(l - 0) [l + 0(1 + 2a)] 


f D (T) 


(A- 11a) 




(l - a) [l + 0(1 + 2a)] 


fj (T) 


(A-llb) 


where we have defined 


f D = 


T 5/2 e' VT Q „1A 2 /»“ A k 

Q Q Q 4 k 

c r v 


3/2 


’ 9 D = h D m M /k (A ‘ l2a) 


= 51, 970 K 


f I = 


t 5/2 - e iL /T 
T e 

Q elA 


3/2 


2TT m E k 


m^/k (A-l2b) 


Equations (A- 11) define a, 0 in terms of p, T, and could be used. However, 
they are difficult to solve for a and 0, and can be simplified considerably 
for hydrogen without much loss of accuracy. Equilibrium hydrogen 
dissociates almost completely before it ionizes significantly. We can thus 
separate the two reactions, and assume a is small during dissociation, 
while 0 goes from 0 to 1. Then 0 is near unity while a increases from 
near zero. These two approximations are separated by a temperature T 


v 


% 
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below which there are molcules and atoms with very few electrons, and 
above which there are atoms, ions and electrons, with few molecules. 

Therefore we use 


* B‘ 

< T s — 


f D ,T) Sa 2 f I (T) 


1 - 


p '1+0 


(A- 13a) 


. - .2 f n (T) 2 fj (T) 

* (1 - a) . D 5 = -= 

T>T : 2d -0)(i + a) " P ' 2 (i - a 2 ) p 


(A- 13b) 


There are easily solved for 0 and a to yield 


T < T*: 0 = (1 + p/f D f /2 > a 


uj i. 

0 p 


1/2 


(A- 14a) 


2 

T > T*: a :: (1 + P/2fj) > 1 - 0 = ^ j + a ^“ 2f D 


(A- 14b) 


In general we will only use 0 as variable below T‘ from (A- 14a), and take 
a = 0. Above T* we will take a as variable, and 0 = 1. The second 
expressions in Eqs. (A- 14) will only be used if an estimate is needed 
for a below T* which differs from zero, or for 0 above T which differs 

from unity. 

To complete the specification of f D and ^ we need the partition 
functions. The combination Q c Q f Q y can be found in Ref. Al, p. 163, 
Eq. (7. 39), where it is called Q yjm . It may be divided as follows: 

Q r = T/ 2 e r . q v 4-.' VT ) (A - 1 
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The factor 2 in has been added to Ref. Al because of the symmetry of 

the hydrogen molecule. The vibrational and rotational temperature constants 

0 and 0 come from Eqs. (7. 32) and (7. 33), p. 162, of Ref. Al using the 
v r 

JANNAF spectroscopic constants (in cgs units) 


B = 60.848, uu = 4405.3, x = 125. 325/4405. 3 -. 0284 
e e e 


6 = 3. 0664/60. 848 = . 0504, y = 60 - 848/4405. 3 = . 0138 


They are 


0 

v 


h p c “e (1 
k 



0 

r 


h P cB e 

k 



0 = 5978 K 

v 


0 = 85. 34 K 

r 


The remaining factor Q c is the coupling partition function. 




3T 0 


. 2x 0 /T 

L + e_v 


r 9 v /T ( 0 / T 

e - 1 l v 




(A- 16) 


e , -2 5. 68x10 2 0 /T 

?n a* . -5 _ 5.04x10 v 

= 1 + 1. 785 x 10 T + — r-yr + ~ 

T V T , / 0 /T 

e - 1 I v 


/ 0 /T V 

l« v -0 


The last expression needed is for the electronic partition function 
of the atoms, which is 


Q elA = 2S 1 6 


9il/t . Sl =t 

j=i 


.2 6 IL /T ^ 

) « 


(A- 17) 
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The first term is the ground state weight 2. The subsequent terms account 

for excited electronic states. The upper limit of the sum must be finite, 

since the terms grow like j 2 for large j. The finite cut-off means that 

above some degree of excitation the electron is no longer considered bound 

to an atom, but is a free electron. Related to this necessity for cutting 

o 

off the sum is the lowering of the ionization potential. The value of h JL 
is not that for removing an electron infinitely far from an atom, because in 
an ionized mixture the electron no longer is considered bound when it is 
far enough from the atom to be indistinguishable from other electrons. 

Thus 0 is not the value 0j = 157, 770 K obtained from the spectroscopic 
constants for a very weakly ionized mixture, but a lower value (hence 
the subscript L). If these two effects are not accounted for, the composition 
of equilibrium hydrogen at high temperatures and pressure will be noticably 

in error, as found in Ref. A2, Table 7. 1. 

The lowering of the ionization potential is discussed in Ref. A3. 

In Eq. (3.85), p. 218, it is suggested that 


3 1/2 

AI = k A 0 = 2 e (2TT n E /kT) 

where e is the charge in esu, and the other units are cgs. For hydrogen, 
e = 4.803 x 10' 10 esu, and we have 

.-20 


AI = 4. 728 x 10 


n E (cm ) 
T(K) 


1/2 


4. 728 x 10 


-30 


n £ (m" 3 ) 

T(K)~ 


1/2 


ergs 


Joules (A- 18) 
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This relates the lowering of the ionization potential to T and the electron 
density, which is in turn related to a by Eq. (A-4), and so depends on 

0 IL throu 8 h f j in Eq. (A- 14b). Thus the determination of a is an iterative 
process, requiring iteration on 0^. It also requires iteration on the 
upper limit L of the sum in Eq. (A-17), which is related to 0^. There 
are a number of such relations which can be used. Here we have chosen 
a simple one, namely to cut off when 

yj 2 > Al/k . 

In the present work we have used one step in the iteration, calcu- 
lating a by the following procedure: 

1) Find a from Eq. (A- 14a) using 0 JL = 6j = 157, 770 K, L = 17. 

2) Find n^, from Eqs. (A-4) and (A-7) using this a. 

3) Find AI from Eq. (A- 18). 

4) Calculate 0 JL = 0j - AI/1.3806 x 10 -22 . 

5) Calculate a cut-off L = greatest integer in - 9 )]^^. 

6) Recalculate a using 0j L and this L. 

The resulting value of a is taken as the correct one. Further iterations 
could be performed, but it will be seen that this one iteration yields quite 
accurate results. 

This calculation of Q #JA is used only for T > T*, where a is the 
varying reaction coordinate and 0 = 1 in general. (However, if an estimate 
of 0 different from unity is needed, we may use this Q in the second of 
Eqs. (A- 14b). ) For T < T , there is very little ionization and the tempe ra- 
ture is low, so we may take Q ^ = c, its ground state value, to calculate 
the varying reaction coordinate 0 . The value of a in this region is zero. 

(If an estimate of a different from zero is needed, we may use the second 
of Eqs. (A- 14a), with Q =2.) 
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The derivatives of a and (3 are needed for the wave and nozzle 
computations. From Eqs. (A- 14) it is easy to find them. In the appro- 
priate temperature ranges we have: 


/ii\ . 8(1 - 8 2 ) dlnf D /3a\ . & Ull 

\ aT /p" 2 dT V T /p 2 dT 

In the differentiation of f , the complicated expression of d In Q /dT is 

D c 

ignored as a small correction, its value taken to be unity. (The value 

of specific heat so derived are shown later to be quite accurate. ) Only 

Q and Q have to be differentiated. The result is 
r v 



( S In z \ 1 / SB \ 

8 In T/p ~ 1+8 \81n T y p 


(A- 19a) 


A similar differentiation of fj involves only the factors T' and S^, since 
the exp (-9j^/T) factor cancels with the one in Q e ^. The result is 
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(A- 19b) 
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z = 2 (l + a) , 




The derivatives with respect to p are even similar, since and fj 
are only functions of T. 



Z = 2(1 + a). 


/ag\ . j (i - g 2 i 

W/ T ' 2 » 

. 6(1 - 8 ) 

In p / T 2 

0 _ a il - a 2 ) 

' \^P/ T " 2p 

/ a In z \ _ q(l - q) 
\ainp/ T " 2 


(A- 20a) 


(A-20b) 



The remaining thermodynamic properties needed are the enthalpy, 
and the specific heat c p . The enthalpy is found by summing the contributions 
of the various components. The translational specific heats of the atoms, 
ions and electrons are 5k/2 per particle. For the ions and electrons, 
this is the total specific heat, so their values per unit mass are 

c pl = 5k/2n V c p£ = 5k/2m E = c pJ n^/n^ (A-2l> 

The corresponding enthalpies are the integrals with respect to T. In 
addition the ions carry the dissociation energy h^° and the ionization energy 
h TT °. Therefore 





The atoms have a translational contribution 5 kT/m^ carry the dissociation 

energy h D ° and also have a contribution from the electronic states derived 

from Q . of Eq. (A- 17) by 
elA 


kT 2 2taQ .lA 


elA m , 


ST 


■*•(-%) 


-where S has been defined in Eq. (A- 19b). Thus the atomic enthalpy rs 

6 


_ 5 kT 
^A 2m, 


+ h D + h IL ^1 - ^ 


(A-23) 


Finally, the molecules have a translational contribution 5 kT/2m^ from 

the two rotational degrees of freedom corresponding to Q = T/20 , and 

r r 

a vibrational contribution found from by 

_ kT 2 * Q v _ kT_ 

v " m M 2T " m M * V 


where ill has been defined in Eq. (A- 19a). Thus the molecular enthalpy 
T v 


is 


kT 


M m 


M 


(i • ..) 


(A- 24) 


The enthalpy of the mixture is the sum of the species enthalpies, 
weighted by their mass fractions: 


* 2 N °i /0 


(A-25) 


The mass fractions are found in terms of a and 3 by using Eqs. (A-4), (A-6) 
and (A-7): 
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(A- 26) 


P Pa Pi Pe gttm E 

P M _ , B _A - mi. ni. — = P a» “T" = 

P P 


-£2-= 1 - P, — =9(1 - a». 

P p 
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I 


use 


I„ th. low temperature region there are only molecule, and atom., we 
h = h M P M /P + h A p A /p with « • #. to h A we ignore the la.t term .mce the 

first t'.rm® of S 2 and S, are equal. We then find from Eq.. (A-23), (A-2 
and (A- 26): 


<T *. h = JJsl. [(0.7 t 0.2* > (X - B) + P (1 + o. 2 e D /T >] (A- 27a) 
- ■ 2m A 


In the high temperature region there are atom., ion. and electron., but 
g = 1. Then Eq.. (A-22). (A-23), (A-25) and (A-26) yield 


* , 5 kT 

T>T : h = 2 

A 


0. 2 0 rt 0.4 6 r s 


(A- 27b) 


The equilibrium specific heat i. obtained by differentiation at constant 
pressure. Two new function, of T enter from the derivative, of * y and 
They are 

6 /T 

= t rn, /dl + i 
V v 


» M 2 e V ■ T d* v /dT + ♦, 

V V 


L 

E 

3 = 1 


S 4 V £ f 2 exp (6 1L /T i 2 . - - (T 2 /8 il ) d S 2 /dT 


(A-28) 


With those, and dS,/dT from Eq. .(A-19b), th. derivative, of (A-27) are 
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T < T*: 




There remains only the specification of T , the temperature at 
which we shift from dissociation with no ionization to ionization with full 
dissociation. Probably the most accurate compilation of the properties 
of equilibrium hydrogen is that of Patch, Ref. A7. A perusal of his tables 
shows that the first temperature at wi.ich the number of hydrogen molecules 
is less than the number of electrons is given by the following table: 

p (atm): 1 3 10 30 100 

T (10 3 K): 7. 3 8 9 10 11 

If we take these values as T , except using 7000 instead of 7300 at 1 atm, 
then the relation 

T* = 2000 log^Q p(a.ti ) + 7000 (A- 30) 

fits p 1, 10, and 100, and is very close to 'he values at p = 3, 30. We 

$ 

will adopt Eq. (A-30) as the definition of T . 
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A comparison of some values calculated using the expressions 
given in this Appendix with those given by Patch is shown in Table A l, 

g 

The enthalpies from Patch have had 2. lb x 10 J/kg added to them because 
of his use of the atom as the zero of enthalpy, rather than He molecule 
which is used here. 

The comparison shows the present properties to be an excellent 
approximation to those resulting from Patch's elaborate calculations. 

The major species are in agreement within a few percent, as are the 

enthalpies and the specific heats. Only the minor species have serious 

* * 

differences, that is, electrons below T and molecules above T . It 
should be noted that Patch has included H^, H^ + , H which are not 
considered here. He also distinguishes the ground electronic state and 

-i- ^ sjf 

excited electronic states of both (H^ and ) and H (H , H ). We 

have added these in our comparison. He also has two cut-offs for Q ; 

6 liV 

one is that used here, the other a distance cut-off. His values of c 

P 

were obtained by numerical differentiation of h. 

Considering the relative simplicity of the present model, it provides 
excellent accuracy for the thermodynamic properties of equilibrium hydrogen. 

Although it is not used in the integration, the Mach number is an 
interesting parameter which can be calculated when p and T are known, 
since they determine the speed of sound ?. For a gas in thermochemical 
equilibrium, the second law of thermodynamics and the expression of h as 
a function of p and p shows that 


2 

a 



(6h/Sp) p 
p" 1 - (dh/*p) 

P 


Conversion of these derivatives to h as a function of T and p, and use of 
(A-l) shows that 
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APPENDIX B 


ABSORPTION OF LASER ENERGY 


The absorption of the incident laser energy by the plasma is depen- 
dent on k L , the absorption coefficient of the plasma. Based on the studies 
of Ref. A5, we use as the absorption mechanism inverse Bremsstrahlung, 
both electron-ion and electron-neutral. The expressions for k have 

J-j 

been given in Ref. A5. The electron-ion absorption is given in Eq. (2. 1) 
there as 

k LEI = °EI n E n I [eXl3 (h P v/kT) ’ 

ct EI is the electron-ion absorption cross-section and v is the laser frequency 
(for 10.6 jim in the present case). The cross- section is given in Eq. (2. 2) 
of Ref. A5 as 


0 


El 


4 / 2tt \ ^ Z 2 e 6 

3 l 3 m kT ) 

\ E / h cm 


1. 63 x 10 


•32 


, 1/2 


(cm 5 ) 


where Ze is the ion charge, e is the electron charge and c is the speed of 
light. Combining these two equations we find 


"LEI 


1. 63 x 10 

m 1/2 


-32 


» E n l( e 1357 ' T . 1Mcm -‘, 


where T is in K and n_, n T are in cm"' 5 . 
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The electron and ion number densities are related to the pressure 
by Eqs. (A-6), (A-7) and (A-8) as 

n = n = 2a 3 p/kT Z (B-l) 

£ 1 

The above expression for ^ is a semi-classical result. To 
account for quantum - mechanical effects we multiply by a Gaunt factor 
of 1. 5, which is an average value for the range of interest. The final 
expression for the electron-ion inverse Bremsstrahlung absorption 
coefficient is then 


‘LEI 


= 5. 1 x 10 


2.2 2 

a P p 

rr.5/2 „2 


, 1357/T 
(e - 1) 


(B-2) 


a = 0 for cgs units, a = 4 for mks units. 

Of course, this expression is only valid where there are ions, a * 0, so 

* 

in our approximation for the chemistry, it is to be used only above T , 
where 8 = 1. 

The absorption coefficient for electron-neutral inverse Bremsstrah- 
lung is also given in Ref. A5, Eq. (2.6), as 


k LEN " °EN n E n N 


[1 - exp (-h p v/kT)] 


where Q is the absorption cross-section and n is the neutral number 
EN w 

density. Equation (2. 8) of Ref. A5 gives 


_ 2 .96 x 10~ 45 T 

°EN " 1 - exp (-h p v/kT) 






(cm 5 ) 


- 80 - 


where { may be approximated by the first two terms of Eq. (2. 9) of Ref. A5, 


when a/T /158, 000 is small, as 


1/2 


{ = 4. 862 (kT/mj hj ) [1 - 0. 2096 (kT/nij hj )]. 


Since for hydrogen h^ /k w 158, 000 K, then we have 


k T _ = 1.6 x 10- 38 JFn„ n exp [-0.0l2zjr (1 - 5. 28 x 10 -4 Jt )] (cm' 1 ) 
LEN E N 


where n^,, n^ are in cm 3 . 

n is given in Eq. (B-l) and the only neutrals present are atoms 
E 

(8 = 1). Their number density is found the same way to be 


n E = 2 a P/ m M * n N = n A = 2 (1 - a) P /kT Z 


Then the final expression for k T E ^ is 


k = 3. 33 x 10 tt (1 - a) B — exp [-0. 0122 Jr (1 - 5. 28 x IQ^*/t )] 
LEN 3/2 _2 1 


b = -6 for cgs units b = -2 for mks units (B-3) 

The expression for k L to be used in the laser absorption expression, 
and in the definition of optical depth T is the sum 


v - v + k 
L LEI LEN 


(B-4) 


A plot of k L£I (dashed lines) and k L£N (solid lines) is given in Fig. B-l 

as k /p vs T. It can be seen that electron-ion absorption dominates at 
L 

temperatures above 11,000 K at 10 atm, above 9, 000 K at 1 atm. The 
sum k L follows the k L£I curves closely above 12, 000 K. 
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appendix c 


THERMAL CONDUCTIVITY OF HYDROGEN 

The only transport property of hydrogen needed is the thermal 
conductivity. Since we deal with hydrogen in thermal equilibrium we 
need no. only the conductivity due to the motion of the particles, usually 
called the frozen conductivity, but also that due to the transport of energy 
of the various species. The total heat flux in a chemically reacting gas. 
in the x direction, can be written as 


- . x & + Y p. u. h. 
" ax Lu i 1 


The first term is the transport of energy relative to the mass average 
speed, and the second term is the transport of energy with the diffusion 
velocity of each species, » , which is the average velocity of the i-th 
species with respect to the mass average velocity. The diffusion velocity 
depends on concentration gradients and temperature gradients. In 
equilibrium, the concentration gradient also depends on the temperature 

gradient, so we may write 


u = - D. dT/dx . 
i i 
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Then the heat flux becomes 


q = -ft F + £ °i D i h i’* T/ * ‘ 


The summation term i.th. contribution caused by the differing average 
velocity of each specie., carrying the enthalpy of that epeci.s, and tncludes 
change, of concentration due to chemical reactions. The total thermal 

conductivity is the sum 


- 83 - 


X c = *F + 2 p i D i h i • 


The calculation of the D. is quite complicated since it involves the 
species diffusion coefficients, as well as the changes in local concentration 
due to reactions. Very little work on this subject has been done for hydrogen. 
The only work known to the authors which covers the high temperature range 
is that of Yos, Ref. A6. Yos provides a table of \ Q for hydrogen at pressures 
of 1, 3, 10 and 30 atm for the temperature range 1000 (1000) 10,000 (2000) 

30, 000. Since duplication of his calculations would be a difficult task, we 
have chosen to use his table and interpolate. The appropriate table (Yos' 
Table III) is reproduced here in MKS units as Table Cl. 

At each pressure, Xq has two peaks, one around 4000 to 5000 K, 
and one around 14000 to 18000 K. The first reflects the effect of dissocia- 
tion, the second the effect of ionization. 

In order to obtain X c for intermediate temperatures, at any of the 
four pressures in the table, a quadratic interpolation scheme was used, 
with the middle value of T taken nearest to the value of T at which X is 
desired. ^ 

For T < 1000, a linear fit was made between 0. 437 at 1000 K and 
0. 177 at 300 K, 


Xc = 3.71 x 10' 4 T + 6.56 j. 10‘ 2 (C-l) 

The value 0. 177 at 300 K was taken from p. 574 of Ref. A7, which gives 
an experimental value of 4227 cal/cm-s-K. 

For the computation of the thermal conductivity loss due to lateral 
temperature variation, the integral 
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TABLE Cl 


EQUILIBRIUM THERMAL CONDUCTIVITY OF 
HYDROGEN X c (W/mK) 

(Yos, Ref. A6, Table III) 


T (10^K) p (atm) 



1 

3 

10 

30 

1 

0.437 

0.437 

0.437 

0.437 

2 

0.903 

0. 862 

0. 837 

0. 823 

3 

6.90 

4.62 

3. 12 

2. 33 

4 

14.8 

15.5 

12.3 

8.77 

5 

3.63 

6.49 

11. 1 

i:. 3 

6 

2.24 

2. 78 

4.40 

7.43 

7 

2. 36 

2.48 

2. v 'J 

4. 00 

8 

2. 71 

2. 71 

2.82 

3. 22 

9 

3. 22 

3. 12 

3. 10 

3. 22 

10 

4.00 

3. 72 

3. 57 

3. 55 

12 

5.90 

5. 31 

4.88 

4.61 

14 

7.58 

7. 11 

6.51 

6. 09 

16 

6.48 

7.67 

7.83 

7. 54 

18 

4.23 

6.20 

7.87 

3.48 

20 

3. 30 

4. 70 

6.81 

8.44 

22 

3.32 

4. 19 

5.85 

7.78 

24 

3.71 

4. 34 

5.60 

7. 37 

26 

4.30 

4. 88 

5.90 

7. 38 

28 

5. 00 

5.57 

6.49 

7. 76 

30 

5.77 

6.38 

7. 33 

8. 52 
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X C (T) dT 



T 



is needed. This was obtained at the same temperatures and pressures 
by a Simpson's rule integration of \ , as interpolated in Table Cl. The 
value at 1000 K is the integral oi (C-l): 




T 

1 ^ ( 1000 ) = J* (:. 71 X lo ’ 4 T + 6. 56 X 10 " 2 ) dT 

300 

= 3. 71 x 10 -4 x . 5 (10 6 - 9 x 10 4 ) + 6. 56 x 10 -2 (1000-300) 
= 215 W/m 


The values so obtained are given in Table C2. At other temperatures, the 
integral was obtained by the same interpolation scheme as used for Xr-* 



TABLE C2 


integral 


OF THERM*. L CONDUCTIVITY 1^ (W/m) 

C 


T (10 K) 


215 

424 

4167 

16606 

25006 

27815 

30096 

32618 

35560 

39148 

49084 

63028 

77280 

87770 

95141 

101699 

108696 

116688 

125976 

136718 


p (atm) 
3 


2734 

2099 

1917 

14452 

10674 

7626 

25005 

22832 

19528 

29356 

30149 

29690 

31942 

33680 

35184 

34522 

36510 

38729 

37421 

39455 

41921 

40825 

42774 

45279 

49820 

51170 

53369 

62447 

62612 

64074 

77565 

77166 

77789 

91440 

93049 

93972 


102175 

107712 

110995 

110955 

120254 

127174 

119420 

131612 

142254 

128615 

143064 

156942 

139045 

155412 

172019 

150947 

169127 

188139 


87/88 


ORIGINAL PAGE IS 
OF POOR QUALITY 




appendix d 


RADIATIVE PROPERTIES OF HYDROGEN 


Radiative transfer influences the hot hydrogen plasma in two ways. 

A portion of the radiation escapes from the plasma and therefore corresponds 
to an energy loss term in the fluid dynamic equations for the system. In 
addition, some of the radiation is transported into he cold gas upstream of 
the LSC .vave and is absorbed there. The rate at which this energy trans- 
port occurs is a primary factor in controlling the propagation of the LSC wave, 
since without the initial heating by either radiation or thermal conduction, 
the hydrogen cannot absorb the laser radiation. In this Appendix we review 

the modeling of the radial radiative losses which was undertaken in the pre- 
A2 

vi ou s report and we introduce a model of the forward axial transport 

Ap 

which is suggested by our previous modeling of LSC waves in air, 

Dl Radiative Loss From the Bulk Plasma 

The radiative losses from the hot hydrogen plasma can be determined 
from a knowledge of the spectral radiative properties of an element of plasma. 
For the temperatures and pressures expected in the hot core of a laser- 
heated rocket engine, that is, temperatures between 1000CK and 20000K and 
pressures between one and ten atmospheres, the plasma is in local thermo- 
dynamic equilibrium. The radiative properties are then summarize- 1 by the 
effective spectral absorption coefficient (T,p) which includes the effect a 
of stimulated emission. 

The calculation of k^ need only include those transitions which make 
significant contributions to the radiation in the regime of Interest, The con- 
centration of hydrogen molecules is negligible in the hot core, hence only 
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transitions involving atoms, ions or electrons are important. Three distinct 
types of transitions occur — free-free transitions (bremsstrahlung), free- 
bound transitions (photo- recombination), and bound-bound transitions (spec- 
tral lines). For our purposes, bremsstrahlung and phott. -recombination into 

the excited states can be adequately described by the semi-classical approxi- 

A9 

mation (Kramer's formula). The quantum mechanical expression for 

Ain 

photo- recombination into the ground state is given by Bethe and Salpeter. 

The strength of an atomic line is usually expressed by the oscillator strength 
f. Quantum mechanical calculations of f are used for the dominant lines, and 
semi- classical formulas for the higher series members. The profile of the 
lines in the Lyman and Balmer series are important for estimating reabsorp- 
tion. The asymptotic Holtt mark profile adequately describes the Lyman 

A1 1 

series lines for this purpose. The exact Stark-broadened profile of Griem 
is required for H , the first member of the Balmer series. Reabsorption 
is unimportant for higher members of the Balmer series. A more detailed 
description of the radiative properties of hydrogen can be found in the pre- 
vious report. Ref. A2. 

A plot of the emission coefficient for a transparent hydrogen 
plasma is illustrated in Fig. Dl. We have chosen a temperature of 18000K 
and a pressure of 3 atmospheres, which are conditions typical of the rocket 
core. The continuum contribution -- bremsstrahlung and photo- recombination -- 
is given by the cross hatched region, whereas the contribution from spectral 
lines is hatched with diagonal lines. For reference, the black body function 
B (T) divided by 2R (for R = 0. 5 cm) is also plotted. The ratio e divided by 

UU _ J uu 

B /2R gives k* (T, p) in cm" . Examination of Fig. Dl reveals that the bulk 
of the radiation occurs in the spectral region of uj greater than 8 eV. This 
radiation is produced by transitions to the ground state. It should also be 
noted that above 13 eV the spectral line contribution no longer shows distinct 
lines but rather resembles the continuum. Therefore the continuum for 
transitions to the bound state effectively extends to photon energies below 
13 eV, a fact which is used in the subsequent modeling of axial transport. 
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B /2R for R 4 0.5 
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Fig. Dl Transparent emission coefficient for 18, 000K, e atm hydrogen. 

Continuum is cross hatched, lines are hatched with diagonal lines. 
Blackbody spectral intensity B , divided by 2R for L = 0. 5 cm, is 
shown for comparison, 111 
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The emission from the Lyman a line swamps all other radiation in 
the transparent approximation. It is essential, therefore, to incorporate 
into the model the effect of reabsorption. The amount of reabsorption de- 
pends upon geometry as well as the pressure and temperature distribution. 

In order to model the bulk radiative loss, the rocket is assigned a charac- 
terise pressure p and a characteristic radius R. (The radiation does not 
change significantly for the small deviations from p and R which occur within 
the rocket core. ) Then the radiation emitted from an element of plasma at 
temperature T can be estimated by calculating the average radiation emitted 
per unit volume from an infinite cylinder of plasma which ha, radius R, tem- 
perature T, and pressure p. The exact solution for the average loss per unit 
volume for any spectral region can be adequately approximated, for the pre. 
sent scenario, by using ,h. optically thin limit for the emission coefficient 
per unit volume whenever the criteria 2^R s 1 i, satisfied, and by using the 
black body spectral function to describe the radiation per uni. surface area 
whenever 2k;R * 1 is true. The portions of the spectrum which are extreme, 
ly black, for example the center of strong lines, should not be considered as 
a bulk loss. Their radiation is controlled by surface effects and cannot be 
reliably predicted without information about the radial temperature gradient. 
Smce these gradient, are not included in the fluid dynamic formulation, w. 
also neglect radiation which is sensitive to these gradients. In the cases 
studied, the spectral absorption coefficient corresponds to the plasma being 
transparent or marginally black WR a 3), or else to the plasma being 
strongly black WR 2 10). In our estimates of bulk radiative losses, the 
strongly reabsorbed contribution (2k' R 4 10) is omitted. 

Numerical calculation, have been made of the bulk plasma loss verm, 
temperature for a variety of pressure, and radii. For insertion into compu- 
ter codes the numerical prediction, can be fitted by an analytic form mg. 
gested from an analysis of the semi-classical continuum radiation. The 
power lost per unit volume P T , is given by 
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P T = 8. 9 x 10 
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where and C 2 are constants, ^ is the ionization potential expressed in 

degrees K, and n^^ is the number density of atoms in the ground state in 

cm . Expressed in terms of the number of atoms, n., we have n 

A AG 

2n A /Q elA’ Where n A is found from (A - 4 >» ( A - 7 ) and ( A -H) of Appendix A, 
and the electronic partition function of the atoms Q is given in (A- 17). 

When expressed in MKS units, as used for calculations, we then find 


17.2 x 10 


n CT 

-9 3 


I («I + ^ Q e 


A_ 

elA 


S 6 I /T 


W 

3 

_m . 


(D-l) 


The values of the constants C 3 and C 2 appropriate to various cases calculated 
ar 9 given in the table in Subsection 2. 3 of the main text. 

D2 Axial Radiative Transport 

An important component of the propagation mechanism for LSC waves 
in air is known to be the forward axial transport of radiant energy from the 
vacuum ultraviolet portion of the spectrum. The same mechanism is expected 
to be important for LSC waves in hydrogen since the hot hydrogen plasma 
emits in the VUV much the same as an air plasma does, and the cool hydro- 
gen atoms upstream of the plasma strongly absorb this radiation. We are 
led, therefore, to model the forward transport in hydrogen in the manner 
which has proven successful for air LSC waves/ 8 

The first step is to determine the region of the spectrum which con- 
tributes to forward transport. For radiation to be an effective transport 
mechanism, it must be strongly absorbed by the cool gas which cannot absorb 
the laser radiation, and it must be emitted strongly by hot plasma. Only 
radiation involving transitions to the ground state atom satisfies both criteria. 
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Fig. Dl can be used to identify the region further. The continuum effectively 
begins somewhere between l2eV and l3eV because of the merging of spectral 
lines. The stronger lines are distinct and should not be included in our esti- 
mate of forward radiative transport since these lines tend to become nar- 
rower at lower temperature, and therefore, even though the lines are strong 
emitters, they do not satisfy the criterion that they also be strong absorbers 
over the same spectral range at low temperature. Morevoer, at lower tem- 
peratures, several of the lines which are shown merged in Fig. Dl become 
distinct, and once again the criteria of strong absorption in the cool gas is no 
longer satisfied. As a conservative estimate of the portion of the spectrum 
which contributes to forward radiative transport, we choose the continuum 
having photon energy cu a 13. OeV. 

An expression for the absorption coefficient for photoionization of the 
ground state has been derived by Bethe and Salpeter (Ref. AlO, pp 304-5), 

In the region of interest in the present calculations (13 < cu < 26eV), they give 

/u) \ 8/3 

\ f cm ^ = 6.34 x 10 18 n AG [cm -3 ] , (D-2) 

where uuj is the photon energy at the beginning of the true continuum, namely 
(Uj = 13. 6eV. The full expression (D-2) can be used in numerical computa- 
tions; however, for discussions of the general properties of the VUV radia- 
tion, it is convenient to define a single absorption coefficient which is charac- 
teristic of the VUV. An appropriate value is 

^VuvIa. ^ 34 x 10 " 18 n AG t cm " 3 ] 

For the computations performed by the computer, the molecular 
hydrogen concentration never gets large enough to warrent its inclusion in 
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the estimate of ky^y However, in order to justify the assertion that all 
radiation emitted in the VUV is absorbed in the cool gas, we have computed 
kyj V for temperatures as low as 2000K. Here absorption by molecular 
hydrogen dominates. The molecular absorption spectrum is characterized 
by an ionization continuum for uu i 15. 5eV, a dissociation continuum for 
15. 5 ^ u) 2: 14. 8eV and sharp absorption bands for uu £ 16. 8eV. While it is 
not easy to associate a single absorption coefficient with the VUV region, 
uu £ 13. OeV, we choose a conservative value of 


( k vuv JB”' 1 ! ■ 


5 x 10" 18 n 


-3-i 


M 


[cm" ] 


where n is the number density of molecule, in order o estimate the effects 
M 

of molecular absorption. This value represents half the experimentally 

Al2 

observed cross section at uu = 17.7eV, as recorded by Hudson. 

A plot of 


k VUV = (k VUV^ A U (k vuv } M 

is shown in Fig. D2 for pressures of one, three, ten and thirty atmospheres. 
In order that all four pressures may be plotted on the same graph, the absorp- 
tion coefficient is divided by the pressure in atmospheres. It can be seen that 
at temperatures expected in the hot core (T ~ 18000K), the mean free path of 
VUV photons is 4 cm at one atmosphere, 0. 6 cm at 3 atm, 0. 1 cm at 10 atm 
and 0.025 cm at 30 atm. Since the hot core does not change temperature 
appreciably over a mean free path in the hot plasma, the radiation stays near 
equilibrium with the plasma, and the radiative transport can be approximated 
by the radiative conduction model that is, the divergence of the VUV flux is 
given by 
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VUV 


dx 
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VUV 


where Sy^jy is the plasma radiative flux in the VUV, and X^y is the radia- 
tive conductivity defined by 


X VUV 


4rr 

3 



VUV 


do) 

k' 

to 


The effective absorption coefficient is found from (D-2) as 


k' 

UJ 


-ou/kT 
(1 - e 



(D-4) 


(D-5) 


where k is the Boltzmann constant. For w* 13. 0 eV and T less than 

30, OOOK, the exponential term is negligible. The derivative of the spectral 

black body function is 


b B 

U) 

5 T 


5.040 (/ e >/T f W 
11,605 t 2 (e (l)/T _ ^2 Lcm 2 eVK 


(D-6) 


where T and oj are measured in eV. 

When (D-5), (D-2) and (D-6) are inserted into the integral of (D-4), 
the integration over u> is quite complicated. To simplify it, we change the 
power in (D-2) from 8/3 to 3, and leave off the factor (1 - e^^) 2 in the 
denominator, of the integral in (D-4), which approximations have counter- 
acting effects. Then the integral is 
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7! e 13 X 11 > 605/T ^ (13 x 11,605/T) X 
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where T is now in K. With this, we can write (D-4) in MKS units as 


X VUV 


0. 775 x 10 30 °elA 


1 -ISO, 865/T 
(11, 605/T) 8 



(150, 865/^* 

A ! 


(D-7) 


This is the form used for in the calculation when the factor given below 
in (D-8) is appended to account for size effects. 

This model suffices for the description of the hot core of a one dimen- 
sional wave. At the wave front, however, the temperature changes abruptly. 
In air LSC waves, the slope of the temperature profile is approximately 2. 5 x 
10^ K/cm at one atmosphere. 3 The radiation conduction model cannot 
give reliable predictions of the local radiative properties in the region of 
the wave since the temperature changes drastically within a photon mean 
free path. However, since the conduction model does conserve energy, it 
can give reasonable predictions of quantities, such as wave velocity and core 
temperature, which depend mainly upon the amount of energy transported 
forward, but do not depend upon the details of the wave front. Therefore, 
the radiative conduction model is still applicable provided that the radiation 
from the hot core is absorbed quickly. This assumption can be justified 
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by an examination of the optical thickness of th« wave front. Assuming 
a linear temperature profile between and T max , the optical depth 

is given by 


x (T ) 
- max 


t(T , T . ) ^ 
max mxn 


/ 


k W T) dx * 


i 

dx f 

dT J 


max 


k vuv (T ^ 


dT 


x (T . ) 
min 


min 


A plot of 


20000K 

T(T) h J* ky uy (T) dT/p 
T 

for four values of pressures p is shown in Dig. D3. The optical depth 
between any two temperatures T min and T max can then be computed from 
Tby evaluating 

p (t(T . ) - T (T ))/ t ' (D_I 

F V v min max V dx 


The slope of the temperature in the wave front is expected to vary 
almost linearly with pressure. Based upon the measurements at atmos- 
phere^ 3 , the maximum slope anticipated in a physical wave is then 

= (2. 5 x 10 4 p[atm])K/cm 
dx 


Consider an LSC wave plasma at three atmospheres and 18000K. If the 
wavefront has the maximum slope, the optical depth between 18000K and 
10000K is found from (D-8) and Fig. D3 to be 0. 72. Between 18000K and 
5000K the optical depth is 2. 1. The fraction of the plasma radiation at 
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Fig. D3 Integral of effective absorption coefficient of hydrogen for VUV 
radiation at 1, 3, 10, 30 atm. Optical depth of LSC wave front ia 
proportional to this quantity. 
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18000K which is absorbed before reaching the place in the wave front where 
the temperature is 10000K can be conservatively estimated as 1 - e°^' = 

0. 5. This is an underestimate because radiation from the hot plasma which 
is not parallel to the axis is absorbed at a faster rate. Similarly we find 
0. 88 of the plasma radiation is absorbed by the point where the temperature 
falls to 5000K. It is reasonable, therefore, to model the forward VUV 
radiation as being totally absorbed and to use the radiative conduction model 
to describe forward radiative transport in the VUV. 

D3 Finite Radius Considerations 

The calculation of bulk losses from the hot plasma includes the effect 
of radius in the reabsorption effects. Only radiation which was transparent 
or marginally black (2k^R £ 3) is included in these estimates. The radiation 
which is strongly black is omitted since it is associated with radial gradients. 
However, for calculations made with a radial conduction loss term, as 
described in Subsection 4. 1, it is possible to include VUV radial losses 
whenever they are not included in the bulk loss term. For the 10 kW and 
5 MW devices, the value of in (D-l) is related to the reabsorption of 
VUV — a value of 2.4 corresponds to the inclusion of VUV, whereas a value 
of 0.4 corresponds to the omission of VUV radial radiation. In the latter 
case, the radiative conduction model can be used to estimate radial losses 
in the VUV by inserting the total thermal conductivity \ c + ^yuv * nto t * ie 
model for radial conduction losses. Of course this modification should 
be made only for the large device — the small device already has the radial 
VUV losses included. 

The finite radius also affects the forward axial transport. If the 
radius is much larger than the mean free path of the VUV photons, it is 
appropriate to use the one dimensional model for forward transport. However, 
if the radius is small compared to the mean free path of the VUV photon, the 
forward axial flux is Beverly limited by geometrical considerations. For 




r 
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of radius 


the axial flux S at the center of the end of 
R, a. .hown in Fig. D4. i« given by 


a 


semi infinite cylinder 


s = tt BjT) [l.2E 3 (*R)] , 


where B^T) i. the apectral black body function, k' 
abaorption coefficient. R i. , h . r adiu., and E 


i« the effective spectral 
exponential integral 


00 
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yx 


dx 


t u that the ,ame reiation,hip ~ - »p«ca. 

depth kR holds for planar slabs of thickness R. ) 

In the spirit of the approximations made in the ^ 4 . 

losses the Kit, u * n the termination of radial 

ransparent approximation is also used to include the 

geometrical effect for forward axial transport; that is, whenever 2k- R * j ia 

true the flux is given by the transparent approximation 2 n B k-R * 

27b 777 7"* Crit6ria ^ 1 *“* body 
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semi-infinite cylinder of radius R. 




In the calculation of the radial losses, 
with respect to temperature 


the integral of \ R = X yuv 



i 


is needed. This was obtained at p = 1,3, 10,30 atm by a Simpson's rule 
integration. The results are shown in Table Dl. At other temperatures, 
the value was found by quadratic interpolation. 


TABLE D1 


INTEGRAL OF RADIATION CONDUCTIVITY 

) 


T(10 3 K) 


p(atm) 



1 

3 

10 

1 

0 

0 

0 

2 

0 

0 

0 

3 

0 

0 

0 

4 

0 

0 

0 

5 

0 

0 

0 

6 

4. 14E-3 

1.39E-3 

4. 3E-4 

7 

2. 12E-1 

7. 07E-2 

2. 14E-2 

8 

4.27 

1.42 

4. 28E-1 

9 

4. 62E1 

1. 53E1 

4. 60 

10 

3. 25E2 

1.07E2 

3. 19E1 

12 

7. 10E3 

2. 23E3 

6.47E2 

14 

8. 69E4 

2.41E4 

6.47E3 

16 

8. 81E5 

1. 94E5 

4. 44E4 

18 

8. 39E6 

1.41E6 

2. 54E5 

20 

6.87E7 

9. 50E6 

1.35E6 

22 

4. 49E8 

5. 63E7 

6. 69E6 

24 

2. 35E9 

2.81E8 

3. 00E7 

26 

1.01E10 

1. 18E9 

1. 20E8 

28 

2. 89E10 

4. 32E9 

4. 23E8 

30 

5.94E10 

1.39E10 

1.33E9 
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(W/m) 


30 

0 

0 

0 

0 

0 

1. 57E-4 
7. 33E-3 
1.44E-1 

1. 54 
1.06E1 

2. 12E2 
2. 03E3 
1.28E4 
6. 42E4 
2. 84E5 
1. 18E6 
4. 58E6 
1. 65E7 
5.47E7 
1. 65E8 
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